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We study the scenario in which a massive particle is thrown into a rapidly rotating Kerr black 
hole in an attempt to spin it up beyond its extremal limit, challenging weak cosmic censorship. We 
work in black-hole perturbation theory, and focus on non-spinning, uncharged particles sent in on 
equatorial orbits. We first identify the complete parameter-space region in which overspinning occurs 
when back-reaction effects from the particle’s self-gravity are ignored. We find, in particular, that 
overspinning can be achieved only with particles sent in from infinity. Gravitational self-force effects 
may prevent overspinning by radiating away a sufficient amount of the particle’s angular momentum 
(“dissipative effect”), and/or by increasing the effective centrifugal repulsion, so that particles with 
suitable parameters never get captured (“conservative effect”). We analyze the full effect of the 
self-force, thereby completing previous studies by Jacobson and Sotiriou (who neglected the self¬ 
force) and by Barausse, Cardoso and Khanna (who considered the dissipative effect on a subset of 
orbits). Our main result is an inequality, involving certain self-force quantities, which describes a 
necessary and sufficient condition for the overspinning scenario to be overruled. This “censorship” 
condition is formulated on a certain one-parameter family of geodesics in an extremal Kerr geometry. 
We find that the censorship condition is insensitive to the dissipative effect (within the first-order 
self-force approximation used here), except for a subset of perfectly fine-tuned orbits, for which 
a separate censorship condition is derived. We do not obtain here the self-force input needed to 
evaluate either of our two conditions, but discuss the prospects for producing the necessary data 
using state-of-the-art numerical codes. 


I. INTRODUCTION 

The cosmic censorship conjecture [T] has over the years 
become a cornerstone of classical general relativity. In its 
weak version it states, in essence, that curvature singular¬ 
ities arising in solutions to the Einstein’s field equations 
must be cloaked behind event horizons, so that they are 
prevented from being in causal contact with distant ob¬ 
servers. Despite being strongly motivated on physical 
grounds, the conjecture’s precise extent of validity re¬ 
mains unclear. A notable counterexample involves finely 
tuned initial conditions [2]. The formulation of the con¬ 
jecture may be refined to exclude such examples [3]. 

In a 1974 paper |4j Wald proposed a simple but pow¬ 
erful framework for testing weak cosmic censorship, us¬ 
ing the gedanken experiment of a particle thrown into a 
Kerr-Newman black hole. If parameters can be chosen 
such that the post-capture mass Mf, charge Qi and spin 
Jf satisfy Mf < + Qf, then a naked singular¬ 

ity would presumably form, in direct violation of weak 
censorship. Whether the equations of classical general 
relativity permit such a process has since been subject of 
much investigation. It is usually assumed that the par¬ 
ticle’s energy and electric charge are much smaller than 
those of the black hole, which then places the problem 
within the realm of black-hole perturbation theory. 

In [3] Wald showed that the over-extremality sce¬ 
nario is ruled out when the configuration is that of a 
pointlike test particle captured by an extremal Kerr- 
Newman black hole. Electrostatic and centrifugal repul¬ 
sion, he showed, would prevent a particle carrying suf¬ 
ficient charge and/or angular momentum from entering 
the black hole. The same conclusion was shown to hold 


true also for a spinning test particle dropped from rest 
at infinity along the symmetry axis of an extremal Kerr 
black hole, with its spin aligned along the axis. In this 
case, it is the repulsion force from spin-spin coupling that 
prevents suitable particles from ever entering the black 
hole. 

However, later work has demonstrated that over- 
extremality is achievable when the initial black hole is 
taken to be nearly extremal—if back-reaction effects on 
the particle’s trajectory are ignored. This was first shown 
by Hubeny [5] for a nearly extremal Reissner-Nordstrom 
black hole, and more recently by Jacobson and Sotiriou 
[B] for a nearly extremal Kerr black hole (“overcharging” 
and “overspinning” scenarios, respectively). The nearly 
extremal Kerr-Newman case was subsequently studied in 
Ref. [7]. In all cases, all orbits identihed as capable of 
driving the black hole beyond the extremal limit lie very 
close, in the relevant parameter space, to the separatrix 
between orbits that are captured by the black hole and 
ones that are scattered off it. In Hubeny’s analysis of 
a radially falling electric charge, electrostatic repulsion 
only marginally fails to prevent the particle from falling 
into the hole: The particle’s radial velocity upon crossing 
(what would have been) the event horizon is proportional 
to the ratio ^ 1 between the particle’s energy and the 
black hole’s mass. The amount of post-capture excess 
charge, Qf — Mf, is found to be quadratic in fj. Similarly, 
in Ref. [B]’s analysis of equatorial-plane captures, over¬ 
spinning particles clear the peak of the effective potential 
barrier with radial velocities cx ?), and the post-capture 
excess spin, Jf — Mf, is quadratic in fj. 

This suggests strongly that back-reaction effects can¬ 
not be ignored and may well change the outcome of the 
gedanken experiment. Heuristically, effects of the (elec- 
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tromagnetic and/or gravitational) self-force enter the 
analysis in two ways. First, the dissipative piece of the 
self-force continually removes some of the particle’s en¬ 
ergy and angular momentum, sending them to infinity 
and down the event horizon in gravitational waves. In 
the Kerr case, dissipative effects may accumulate as the 
particle “hovers” above the peak of the effective potential 
on a nearly circular orbit. Second, the conservative piece 
of the self-force might supply just the right amount of 
additional repulsive force to prevent would-be overcharg¬ 
ing/overspinning particles from ever entering the black 
hole. For particles sent in from infinity in the Kerr case, 
this second effect may be formulated in terms of a shift in 
the critical impact parameter for capture: If the gravita¬ 
tional self-force (GSF) shifts the critical impact parame¬ 
ter inward by a sufficient amount (for a given energy-at- 
infinity), then would-be overspinning particles may end 
up being scattered away rather than captured. 

There have been several recent attempts to quan¬ 
tify the effect of back-reaction in the problem. Focus¬ 
ing on the Reissner-Nordstrdm case, Isoyama, Sago and 
Tanaka [5] argued that the full effect can be properly 
taken into account by considering the quasi-equilibrium 
configuration of a charged particle placed precisely on 
the capture-scatter separatrix. An exact solution is 
known for this configuration—the static double Reissner- 
Nordstrom spacetime—and the authors calculated that 
its total energy is always greater than its total charge. 
They have also established that radiative losses during 
the final plunge are negligible, hence concluding that (un¬ 
der the assumption that the true capture system does 
indeed go through a quasi-equilibrium state) the final 
configuration cannot be a naked singularity. 

In a later work, Zimmerman, Vega and Poisson [5] took 
up the challenge of directly calculating the charged par¬ 
ticle’s trajectory including the full effect of the electro¬ 
magnetic self-force. Analyzing numerically a large sam¬ 
ple of orbits within the domain identified by Hubeny, 
the authors found no example of successful overcharg¬ 
ing: All particles with a combination of charge and en¬ 
ergy suitable for overcharging the black hole were found 
to be repelled before reaching the horizon. This anal¬ 
ysis, however, neglected the potentially important ef¬ 
fect of back-reaction from the gravitational perturba¬ 
tion sourced by the particle’s electromagnetic energy- 
momentum. A complete analysis would require calcu¬ 
lation of the corresponding GSF, but techniques for cal¬ 
culating self-forces in the coupled problem are only now 
starting to be developed mm- 

In that respect, the Kerr setup provides a cleaner en¬ 
vironment, in which the perturbative problem is purely 
gravitational (at the obvious cost of abandoning spheri¬ 
cal symmetry). Barausse, Gardoso and Khanna [TH [13] 
studied the dissipative GSF effect in the Kerr over¬ 
spinning problem, focussing on ultra-relativistic particles 
on equatorial orbits. Using analytic arguments backed 
by a numerical calculation of the energy and angular- 
momentum carried away in gravitational waves, they 


showed that dissipation averts the overspinning for some 
but not all of Jacobson-Sotiriou’s orbits. For sufficiently 
small fj, the dissipative effect is always negligible and can¬ 
not prevent overspinning. This result highlights the im¬ 
portance of accounting for the full effect of GSF. To reach 
a definitive conclusion necessitates an actual calculation 
of the full local GSF acting on the captured particles. 

In the past few years, rigorous methods for GSF calcu¬ 
lations in Kerr spacetime have advanced enough to allow 
a more systematic and complete treatment of the over¬ 
spinning problem. The program initiated with this pa¬ 
per revisits the problem from this new vantage point. It 
seeks to obtain a more conclusive answer to the ques¬ 
tion of whether it is indeed the self-force that provides 
the mechanism by which black holes protect themselves 
from being overspun. 

Our current paper lays the necessary groundwork. 
Goncentrating on equatorial orbits, we first identify the 
complete “window” in the parameter space in which over¬ 
spinning occurs if the GSF is ignored. We then formulate 
a condition for this window to be eliminated by the ef¬ 
fect of the full GSF. The condition takes the form of 
an inequality that is required to hold for each member 
of a certain 1-parameter family of geodesics, and it in¬ 
volves the GSF calculated along such orbits. Here we 
do not obtain the necessary GSF data, but we discuss 
methods for computing it (numerically) using existing 
codes. With collaborators we have began work to obtain 
the GSF data, and we intend to present the results in a 
follow-up paper. 

The rest of this introduction summarizes our analysis 
(also in relation to previous work) and describes its main 
results. 


A. This work: overview and results 

Jacobson and Sotiriou [5] assumed that overspinning 
occurs if two conditions are met: (i) the geodesic trajec¬ 
tory of the test particle is timelike at the horizon, and 
(ii) Jf > Mf. The first condition is very lax. It allows 
for low-energy orbits that are deeply bound to the black 
hole and confined to the immediate neighborhood of the 
horizon. The physics of such orbits becomes very sub¬ 
tle, especially when self-gravity and finite-size effects are 
included [consider that deeply-bound orbits below the 
innermost stable circular orbit (ISGO) plunge into the 
black hole within an amount of proper time that shrinks 
to zero in the extremal limit Ml so one would prefer¬ 
ably avoid such orbits as candidates (see, however. Hod 
m for a heuristic treatment). Jacobson and Sotiriou 
acknowledge this issue, and to address it they supple¬ 
ment their analysis with two specific numerical examples 
of overspinning orbits that are sent in from afar. They 
stop short of determining the full range of overspinning 
orbits when deeply bound orbits are disallowed. 

Our capture condition will be more stringent, and more 
in the spirit of Ref. [3] : We will send our particle in from 
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“sufficiently far” (this condition will be made precise in 
the next section), and deem it “captured” if it has no in¬ 
ner radial turning point outside the black hole. Thus, for 
a legitimate capture we demand that the particle “clears” 
the peak of the effective potential on its inward journey. 

In Sec. m] we revisit the overspinning problem in the 
geodesic approximation. We identify the precise region in 
the parameter space of equatorial orbits around a nearly- 
extremal Kerr black hole (excluding deeply bound orbits) 
where overspinning occurs, and give analytic expressions 
for the boundaries of that region. The overspinning win¬ 
dow is illustrated in Figs. and Perhaps unexpect¬ 
edly, we find that only particles sent in from infinity are 
capable of overspinning the black hole. This fact has 
somehow gone unnoticed in previous work, to the best of 
our knowledge. We find that for any given value of the 
particle’s energy at infinity, there exists an open range 
of orbital angular momenta and particle’s rest masses for 
which overspinning occurs. That only orbits coming from 
infinity are potential overspinniners is somewhat fortu¬ 
itous, because for such orbits it is straightforward to iden¬ 
tify the system’s total [Arnowitt-Deser-Misner (ADM)] 
energy and angular momentum even when GSF effects 
are included. 


We then turn to analyze the GSF effect. In Sec. m we 
first review essential results from GSF theory, and then 
discuss the determination of the “critical orbit” that sep¬ 
arates (in the relevant parameter space) between plung¬ 
ing and scattered orbits. We do this in two steps. First, 
we ignore the dissipative piece of the GSF, and calculate 
the correction due to the conservative GSF to the criti¬ 
cal value of the angular momentum for a fixed value of 
the energy-at-infinity. Then we restore dissipation and 
consider its consequences. Under the full GSF, all criti¬ 
cal orbits merge into a “global attractor” that takes the 
system adiabatically along a sequence of quasi-circular 
unstable orbits ending at the ISCO, before plunging into 
the black hole. By fine-tuning the initial value of the 
angular momentum (for a given initial energy), an orbit 
can be made to evolve arbitrary far along the global at¬ 
tractor. We make a formal distinction between “generic” 
and “fine-tuned” captured orbits, based on how the dif¬ 
ference between the initial and final values of the par¬ 
ticle’s specific energy scales with the particle’s mass 
(in a procedure whereby fi is taken to zero while holding 
fixed the initial specific energy and angular momentum). 
Generic orbits are ones for which that difference vanishes 
for p, —>■ 0 (this includes, for example, all of the orbits 
considered in Refs. min]) ; fine-tuned orbits are ones 
for which the difference does not vanish even for /i —>■ 0. 


With this preparatory work in place, we move on, in 
Sec. |IV[ to obtain the overspinning condition as modified 
by the full GSF. The end result are two inequalities, one 
for generic orbits [Eq. ( [^ ] and another for fine-tuned 
ones [Eq. (81)], which describe conditions for overspin¬ 
ning to be averted under the effect of the full GSE. In the 
generic case, the condition involves only the conservative 
piece of the GSE, evaluated along critical geodesics in 


the extremal Kerr limit. Overspinning can be ruled out 
if and only if the condition is met for each member of 
this one-parameter family. The condition for fine-tuned 
orbits requires, in addition to the conservative GSE, also 
knowledge of the fluxes of energy and angular momen¬ 
tum radiated to infinity by particles on unstable circular 
orbits, in the extremal Kerr limit. Overspinning can be 
ruled out if and only if the condition is met for any values 
of the initial and final energies. 

In Sec. we propose an alternative form of the over¬ 
spinning conditions, based on the framework of the “first 
law of binary black-hole mechanics”, as recently applied 
to orbits in Kerr m- The alternative form, given (for 
the generic case) in Eq. (104), involves only perturba¬ 
tive quantities calculated along (unstable) circular or¬ 
bits, which should be more easily computable with ex¬ 
isting GSE codes. This simplified formulation relies on 
explicit expressions, given in EH- for the ADM-like en¬ 
ergy and angular momentum of circular-orbit configu¬ 
rations in Kerr, including leading-order self-interaction 
terms. The underlying theoretical framework is yet to be 
firmly established a nd te sted, however, so we regard the 
simplified condition (104) as somewhat less rigorous than 
our direct condition ( |69p . Our stance is that it would be 
desirable to evaluate both forms of the condition, for the 
sake of establishing confidence in the result. 

In Sec. |VI| we discuss the numerical input required for 
evaluating our overspinning conditions, and the prospects 
for obtaining it through adaptation of existing codes. 
Evaluation of the direct conditions (69) and (81) in¬ 
volves GSE calculations along unbound orbits on a 
nearly-extremal Kerr background, which has not been 
attempted so far. However, we think the basic computa¬ 
tional infrastructure for such calculations is well in place. 

Section [VII| summarizes our results and speculates on 
what a numerical evaluation of our censorship conditions 
might yield. The Appendices contain some of the details 
of calculations done in Secs. nn and lYl a derivation of 
the ADM energy and angular momentum for the system 
under consideration; a calculation of the GSF-induced 
shift in the critical value of the angular momentum; and 
an evaluation of radiative loses during the final plunge 
into the black hole. 

Throughout this paper we set G = c = 1 and use the 
metric signature (—, -|-, -|-, -f). 


II. OVERSPINNING ORBITS IN THE 
GEODESIC APPROXIMATION 

Our initial configuration features a Kerr black hole of 
mass M and angular momentum J = a M < A 

pointlike test particle of rest mass /i ^ M is sent in on 
a geodesic of the background Kerr geometry. As in [5], 
we restrict attention to prograde orbits in the equatorial 
plane, so that the orbital angular momentum is aligned 
with the spin of the black hole. (Intuitively, this con¬ 
figuration seems most favourable for a successful over- 
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spinning.) We denote the particle’s specific energy and 
angular momentum by E and L, respectively; these are 
constants of the geodesic motion. For the geodesic ap¬ 
proximation to make sense, we must assume pE <?; M 
and piL J. Then, clearly, overspinning could only be 
possible, in principle, if the black hole is nearly extremal. 
We write 


a/M = 1 - e^ 


( 1 ) 


Substituting from Eq. Q and solving for E and L in 
terms of the circular-orbit radius, r = R, gives E = 
EfiR) and L = LfiR), with 


E,{R) 

La{R) 


1 - 2R-^ + h.R-3/2 
Vl - -f 2h7?-3/2 ’ 
.Ri/2(l-2h.R-3/2 + a2.R-2) 
y'l - 3.R-1 -f 2h.R-3/2 


( 6 ) 

(7) 


where e ^ iQ 

Below we study the overspinning scenario in the above 
setup, but we begin with a survey of some essential prop¬ 
erties of timelike equatorial geodesics of the Kerr metric. 


A. Relevant results for Kerr geodesics 


Let denote the particle’s four-velocity. In Boyer- 
Lindquist coordinates {t, r, 9, (j)} we have = 0, and 

ut = 0 , = 0 , ( 2 ) 

where an overdot denotes differentiation with respect to 
proper time. The two equalities express the conservation 
of energy E = = —Ut and angular momentum 

7- = where C(“) := df and ^ := d'^ are 

Killing vectors associated with the time-translation and 
rotational symmetries of the Kerr background. The pair 
{E, L\ parametrizes the family of equatorial geodesics 
(up to initial conditions). 

The normalization UaU°‘ = — 1 now gives the radial 
equation of motion, which we write in the form 

= B{r) [E - V. (L, r)] [E - K+(L, r)]. (3) 

Here r is the Boyer-Lindquist radius of the orbit, B{r) := 
1 -I- a^(r -I- 2M)/r^, and (for MaL 0) 


V±{L,r) 


2MaL ( 1^ ^ Br^[L^{r-2M) + rA] 

Br^ I 1 ± V ^ 


(4) 

with A := — 2Mr + a?. For prograde orbits, the po¬ 

tential V- is manifestly negative definite, so the factor 
B{r){E — V-) in Eq. ([^ is manifestly positive definite. 
Thus, plays the role of an effective potential for the 
radial motion, which is allowed for E > V+{L,r), with 
an equality identifying radial turning points. 

Stationary points of V+{L,r) outside the black hole, 
when they exist, correspond to circular orbits. These 
satisfy the simultaneous conditions 


E = V+, = 0 (circular orbits). (5) 


^ Note, to avoid confusion, that has instead a/M = 1 — 2e^. 


Here an overtilde denotes a-dimensionalization using M, 
i.e., R R/M, a := a/M and L := L/M] we shall adopt 
this notation throughout the rest of the paper. Timelike 
circular orbits exist only for R > i?ph(a), the radius of 
a photon’s unstable circular orbit (“light ring”). i?ph(a) 
is the (unique) root of 1 — -f 2aR~^/‘^ greater than 

the event horizon’s radius, Rc\i{a) = 1 -h (1 — The 

angular velocity H := u'^/u* of any circular geodesic orbit 
reads 


n{R) = {d + ( 8 ) 

The number of stationary points of and their loca¬ 
tion depend on L. There are none outside the black hole 
when L is below a certain critical value Lisco(a), and there 
are two for L > Lisco(n): a maximum representing an un¬ 
stable circular orbit, and, further out, a minimum rep¬ 
resenting a stable one. The critical value Lisco(a) marks 
the innermost stable circular orbit (ISCO). It is given by 
Lisco = Lc(i?isco), where the ISCO radius i?isco is found 
by solving Eqs. ^ simultaneously with d‘/.V+{L,r) = 0. 
The ISCO may also be said to represent the outer bound¬ 
ary of the region of unstable circular orbits. 

The radii of unstable circular geodesic orbits span the 
interval i?ph(a) < R < i?isco(«)- This I-parameter family 
of orbits will feature dominantly in our analysis, because 
it defines the capture-scatter threshold where much of 
the relevant physics occurs. Members of the family may 
be parametrized by either E or L, both being mono- 
tonically decreasing functions of R between i?ph (where 
E,L ^ oo) and i?isco for any a < 1. This monotonicity 
can be readily established from Eqs. and Q. Hence, 
the radius R itself is also a valid parameter. 

To each unstable circular orbit there correspond non¬ 
circular homoclinic-like geodesic orbits m that join the 
circular orbit asymptotically in either their infinite past 
or their infinite future, or both. Nearly-homoclinic or¬ 
bits exhibit a “zoom-whirl” behavior [T5]: an episode 
of prolonged rotation (“whirl”) about the location of 
the associated unstable circular orbit. We will see that 
all orbits relevant to the overspinning problem fall in 
that category. Based on the correspondence with homo¬ 
clinic orbits, unstable circular orbits may be divided into 
“bound” {E < 1) and “unbound” {E > I). The radius of 
the innermost bound circular orbit (IBCO) is obtained 
by solving EfiR) = I, giving .Ribco = [1 -f (1 - 

Figure illustrates the range of stable and unstable 
circular orbits, and the location of the various special 
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FIG. 1. Timelike circular equatorial geodesics around a 
nearly extremal Kerr black hole, shown here for a = 0.99M. 
The plot shows specific angular momentum versus Boyer- 
Lindquist radius. Orbits with r > Risco (magenta) are stable, 
while these with r < Risco (blue) are unstable. Also indicated 
are the innermost bound circular orbit (IBCO, E = 1) and 
the photon orbit (“light ring”, E,L ^ oo). In the extremal 
limit, a ^ M, the radii 7?isco, I^ibco and Rph all coincide with 
the horizon radius Rah- 


orbits mentioned, in a particular example (d = 0.99). 
We note the ordering 

Rch ^ .^ph ^ .^ib CO ^isC0 5 (9) 

which applies for any d < 1. 

Let us now specialize to a near-extremal Kerr back¬ 


ground with spin as in Eq. ([^. One finds 

i?oh = 1 + 's/s e + 0{e ^), (10) 

I?ph = l + ^/^e + 0(e2), (11) 

^ibco = (1 + e-f (exact), (12) 

= 1 + + 0(e4/3). (13) 

The function E^^R) in Eq. (§ can be inverted perturba- 
tively in e to obtain the radius of an arbitrary unstable 
circular orbit in terms of its energy E. We find 

R=l + epi{E) + e^P2{,E) + 0{^), (14) 


where the first two coefficients, needed below, read 
2y/2E 2{2E^ -E^ + 1) 


Pi = 


V3E2 _ 1 ’ 


P2 = 


(3E2 - 1)2 


(15) 


Equation (121 is the special case of (14) with E = 1, 


giving Pi = 2 and p 2 = 1. 

It follows that, in the extremal limit e —?► 0, the Boyer- 
Lindquist radii of the light-ring and the ISCO both co¬ 
incide with the horizon radius, and so do the radii of all 


unstable circular orbits enclosed between them. Also pe¬ 
culiar is the fact that the ratio of coordinate differences 
{Risco - Rch)l{Rihco - Rch) diverges as e 0. A closer 
look reveals [T^ that the light ring, IBCO and ISCO re¬ 
main separated from the horizon, and from each other, 
when examined on a Boyer-Lindquist <=const slice: On 
that slice, the proper radial distance between the light¬ 
ring and the horizon is finite, and so is the distance be¬ 
tween any fixed-E unstable circular orbit and the light 
ring. The proper radial distance between the ISCO 
and any fixed-E unstable circular orbit diverges on the 
t=const slice; the geometry of the t=const hypersurface 
appears to “stretch” infinitely around the ISCO location 
m- The situation, however, is rather different when ex¬ 
amined on a horizon-crossing time slice. As emphasized 
recently by Jacobson mi, on any such slice, the light 
ring, IBCO and ISCO all actually coincide with horizon 
generators. From that perspective, they—and all unsta¬ 
ble orbits in between them—are “at the same place” in 
the extremal limit. 

These subtleties will not affect our analysis directly: e 
will be kept small but nonzero, and the strict ordering 
© will therefore apply on any time slice. However, we 
must take note of the degeneracy of i? as a parameter for 
unstable circular orbits when e —>■ 0. The energy E, on 
the other hand, remains a good parameter even in this 
limit, spanning the entire range oo > E > ^. We will 
thus generally adopt E for labelling unstable circular or¬ 
bits. Given E, the angular momentum La{R{E)), which 
we henceforth write as Lc{E), is obtained by substitut¬ 
ing Eqs. ([^ and (141 in Eq. ([^ and then expanding in e. 
The result is 


LaiE) = 2E+{6E^-2)^/^e + Oie^). (16) 

We note that to determine the 0{e) term here required 
the explicit values of both pi and p 2 of Eq. ([T^. 


B. Exclusion of deeply bound orbits 

Heuristically, if we assume our point particle represents 
a compact object—say, a Schwarzschild black hole—then 
its effective proper “diameter” is ~ p,. Below it will be¬ 
come clear that a successful overspinning requires /r ^ e, 
and so relevant objects have proper diameters ~ e. Now 
consider placing such an object in a deeply bound or¬ 
bit with an outer turning point at r < i?isco [and with 
L > Lc{E)]. Such an object (it can be checked) will 
plunge through the horizon within a proper time of 0{e) 
(at most), comparable to its own “light-crossing time”. 
It is not clear whether the object can be made to initially 
“fit” in its entirety outside the hole. At the very least, it 
is not clear if the simple model of a point particle and a 
stationary horizon provides a faithful description of the 
physics in this case. 

To avoid such subtleties, we wish to exclude deeply 
bound orbits from our analysis. We achieve this by re¬ 
quiring that, if the orbit possesses an outer radial turning 
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point at some r = Tout, then 

^out ^ -^isco(^ )■ (17) 

It can be checked that, under this condition, the proper¬ 
time interval along any timelike equatorial geodesic con¬ 
necting r = Tout to r = is finite (nonzero) even in 
the limit e —>■ 0 (taken with fixed E,L). The condition 
0 demands that eligible particles must clear the peak 
of the effective potential (when such a peak exists) as 
they plunge into the black hole. 


C. Overspinning domain 


Given the restriction (17), a necessary and sufficient 
condition for a falling particle of specific energy E to be 
captured by the black hole is 


L < L,{E). 


(18) 


A captured particle would overspin the black hole if and 
only if 


{M + nEf < aM + nL. (19) 

Using a = 1 — and introducing the small mass ratio 
77 := ^/M, this condition becomes 

+ r]W + r]^E^ < 0 , (20) 

where we have introducec 0 

W:=2E-L. ( 21 ) 


Note that Eq. (20) sets a lower bound on L (for given 
E,ri,e), while Eq. (18) sets an upper bound. Also note 
that Eq. (20) implies the necessary condition W < 0 for 
overspinning to occur. 

Our goal now is to identify the complete domain in 
the space of {r],E,L} for which the conditions (18) and 
( 20 ) are simultaneously satisfied, assuming e <C 1 and the 


condition ( [T7| . For easy reference, let us call this domain 
“OS”, for “overspinning”. 

We first show that orbits with L < Lisco all fall outside 
OS. To this end, consider first the ISCO itself, where 
W = 2Aisco-Asco =: Wisco- Using Eqs. ([^ and (1^) 
we obtain Wjsco = —ce‘^^^ + 0{e‘^), where c = 2^/^-\/3 > 0 . 
Thus, lUisco is negative as required, but it can be easily 
checked that ( 20 ) is always violated for sufficiently small 
e: Replacing W —>■ —in Eq. (20) and considering 
the left-hand side as a quadratic function of 77 , we find 
this function is positive definite for any e < (2Eisco/c)^. 
[Since Eisco is bounded from below by Eisco(e ~ 0) “ 


^ Heuristically, iy/2 may be interpreted as the specific energy in 
a co-rotating frame with Q = 1/2, i.e., the common angular 
velocity of all unstable circular geodesics in the extremal limit. 


we find that (20) is always violated for e < ^.] This rules 
out the ISCO itself, and it immediately rules out also all 
orbits with {E > E[sco, L = Tisco}) for which W > Wisco- 
Orbits with {E < Eisco, T^isco} can potentially satisfy Eq. 
( 20 ), but they are always deeply bound in the sense of 
failing to satisfy Eq. ( 0 : For any E < Ejsco, the orbit 
has an outer radial turning point at a radius rout < 7?isco- 

The upshot is that orbits with L = Lisco are all outside 
OS. For orbits with L < Lisco we would need to require 
E < Eisco in order for W to be sufficiently negative. 
But, once again, such orbits are excluded on account of 
their being deeply bound. We conclude that orbits with 
L < Eisco are all outside OS. 

Let us focus therefore on orbits with L > Eisco- For 
such an orbit to be in OS, we require that (given E, 77 , e) 
E is bounded from above by Lc{E) and simultaneously 
from below via Eq. (20): 


-I- 2r]E - 


-77^E^ <riL < r]Lc{E-,e). 


( 22 ) 


We have made here the e dependence of Ec explicit, for 
clarity. The span of the permissible range is 77 A r. := 
—e^ — r][2E — Lc[E\ e)] — rf'E'^, or, using Eq. (16), 


77Ai = —-I- ■qe's /— 2 — rj'^E^, 


(23) 


where we have omitted terms of 0{r]€^). OS is populated 
if and only if we can find E, 77 , e for which > 0. 

A few conclusions can be drawn immediately. First, 
considering tjAl in Eq. (23) as a quadratic function of 77 , 


we find it has a maximum value 


max 77 A/, = 

V 


^(E2 - 1) 

2E2 ■ 


(24) 


This is positive only for E > 1. Therefore, all orbits with 
E < 1 fall outside OS. Bound orbits cannot overspin. 

Second, for any E > 1, we can obtain A^ > 0 by 
choosing the mass ratio 77 from within the interval 


e 77 _(E) < 77 < e 77 +(E), (25) 


where 


V± = 


72 E2 L 


73E2 _ 1 ± 7 t ;2 _ ij . (26) 


In other words, overspinning can he achieved for any E > 
1, as long as 77 satisfies (25). Since the condition Aj;^ > 0 


is both necessary and sufficient, the converse also holds: 
All orbits in OS satisfy E > 1 with Eq. (25). 

Third, from Eq. 


it follows that 77 must be chosen 
to be of 0(e) (assuming E ^ 1/e). One can check that 
77 + has a maximal value of 


max 77 + = VW 2 , 
E 


(27) 


obtained for E = 2/73. Therefore, the range 77 > 7^/2 e 
lies outside OS. The bandwidth of admissible mass ratios, 
for given E and e, is 

A^ := 677+ - e77_ = e72(E2 - l)/E^, (28) 
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FIG. 2. Domain of mass ratios r] for which overspinning is 
possible in the geodesic approximation. 77 is shown divided 
by the near-extremality parameter e = (1 — and E 

is the particle’s specific energy. The boundaries r]±{E) are 
given in Eq. (26 1 . Overspinning is not possible for i? < 1 or 
77 > y/ 3/2 e. However, for any value E > 1 there is a range 
of rj for which the black hole may be overspun. This happens 
if the particle’s angular momentum is chosen from within the 
range indicated in Eq. (291. 


FIG. 3. The overspinning window, shown in the plane of 
E, W (where W = 2E — L/M) for several values of 77 /e. Note 
W is shown divided by e. The boundaries W± are given in Eq. 
(301. The lower boundary W-(E) (which does not depend on 
77 ) arises from the requirement that the particle is captured by 
the black hole. The upper boundary IF+(i?, 77 /e) comes from 
the requirement that the final object is an over-extremal black 
hole. Overspinning is possible with any E > 1, provided 77 is 
chosen from within the range shown in Eq. (251. 


which is maximal for B = \/2. Figure depicts the 
permissible range of rjje ss & function of E. 

Fourth, from Eqs. (l24| and (28) we learn that an E = 


const(> 1) slice of OS has maximal dimensions 
e ^/77 ^ e and ~ e. OS is thus a narrow “tube” in 
the {E, L, 77 } parameter space, of a cross section ^ e x 
e, whose boundary is tangent to the surface of unstable 
circular orbits, L = Lc{E). 

To summarize, we have found that OS is a narrow tube¬ 
like region of the {E,L,t]} space, described by > 1, 
Lc{E-,e) - AL{E,ri;€) < L < Lc{E;e) and er]_{E) < rj < 
£ 77 + (if), where A + and 77 + are given in Eqs. (23) and (26 1 , 
respectively. A neater description of the OS window is 
obtained in terms of the quantity W defined in Eq. (211: 
Rearranging Eq. (22) and using (16), we find 


III. SELF-FORCE PRELIMINARIES 

Because the width of the overspinning window is of 
0(77), self-gravity effects may potentially close this win¬ 
dow, and they must therefore be included in the analysis. 
Specifically, the GSF modifies the capture condition (18) 
by changing the functional relation Lc{E) at 0(77). It 


also modifies the overspinning condition (19) by dissipat¬ 


ing away some of the system’s initial energy and angular 
momentum. In this section we introduce relevant results 
from the theory of self-forced motion. In Sec. |IV| we will 
then use these results to derive conditions for capture 
and overspinning under the full GSF effect. 


A. Equation of motion with self-force 


eW-iE) < IF < eIF+(F;, 77 /e), 


(29) 


where 


IF_ = -\/6E2 -2, 


W+ = -l- + ^E 


(30) 


This domain is illustrated in Figure]^ for a sample of rj/e 
values. To overspin a black hole of given M and e ^ 1, 
one should pick an E greater than 1, choose any 77 from 
within the interval (25), and then choose IF (hence L) 


from within the interval (291. 


There now exists a rigorous formulation of the equa¬ 
tions of motion for compact objects in curved spacetime, 
valid through first post-geodesic order in perturbation 
theory—see and references therein, and [531 HI] 

for recent reviews. The formulation applies in situations 
where all lengthscales associated with the compact ob¬ 
ject are much smaller than the typical curvature radius 
of the background geometry. The motion of the compact 
object is then determined via a systematic procedure of 
matched asymptotic expansions, and interpreted as an 
accelerated motion in the background spacetime, subject 
to an effective GSF (oc 77 ^). One of the results is that the 
internal structure of the object does not affect the self- 









































acceleration at 0 (? 7 ) (except, if the particle is spinning, 
through the familiar Mathisson-Papapetrou spin term). 

The GSF formalism should be applicable in our setup, 
since we work under the assumption <?; M. The 
introduction of the small background-related parameter 
e should not pose a problem, because the background’s 
curvature radius remains much larger than fiE even in 
the limit e —> 0, and even as the particle approaches the 
horizon. We will indeed proceed under the assumption 
that the standard first-order GSF formalism is applicable 
anywhere along the particle’s trajectory until it crosses 
the horizon. 

The equation of motion, including the leading-order 
GSF, may be written in the form 

= F°‘. (31) 

Here u°‘ is the particle’s four-velocity, tangent to the (ac¬ 
celerated) trajectory in the background spacetime (Kerr, 
in our case) and normalized using 

= -1, (32) 


where gap is the background (Kerr) metric. The covari¬ 
ant derivative in (31) is taken with respect to gap, and 
is the first-order GSF, proportional to The GSF 
is normal to the four-velocity, gapu'^F^ = 0, so that the 
rest mass g remains constant. Methods to compute in 
Kerr spacetime are reviewed in |25) . In should be noted 
that F“ itself is a gauge-dependent notion: A full, gauge- 
invariant information about the motion is contained only 
in the combination of the GSF and the metric perturba¬ 
tion with which it is associated [26]. 

Now consider a particle sent in along the equator of 
the Kerr black hole, i.e. with 0 = 7r/2 and u® = 0 at the 
initial moment, where hereafter r is proper time along 
the self-accelerated orbit in gap- In any reasonable gauge, 
the component would vanish from symmetry and the 
motion will remain equatorial. Let us then define 


E := -lit, 


L .— 


(33) 


in analogy with E and L of the geodesic case. Here, how¬ 
ever, E and L are not constants of the motion. Rather, 


Eq. (31) tells us they evolve (slowly) according to 


dE 


^dT " 


(34) 


where Fa = gapF^. With these definitions, Eq. (32) 
produces the radial equation of motion 


f2=i?(r)(K-K_(L,r)) (^E - V+{L,r)) , (35) 

whose form is identical to that of Eq. ([^ —except that 
here E and L are slow functions of r along the orbit. 

The results of the previous section lead us to focus 
attention on particles sent in from infinity, i.e., ones with 
r(T —>■ —oo) —>■ oo. For such particles, we define 


Eoo := E(r-oo). Loo := T(t- t-oo). (36) 


From Eq. (34) we have 


E{j)=Eao+/^E{T), 


where 

fj-AE{T) 



L(r) = Loo+AL(T), (37) 
/iAL(T) = f Fp dr. (38) 

^ —OO 


In principle, the coupled set (35) with (37) determines 
the self-accelerated orbit, given the initial values Loo, Loo 
and a method for calculating the GSF along the orbit. 


B. Dissipative and Conservative pieces of the 
self-force 


The quantities AE{t) and AL(t) encapsulate both 
conservative and dissipative effects of the GSF. This ter¬ 
minology refers to a splitting of the GSF in the form 

L“ = (39) 


where the first and second terms are the self-forces ex¬ 
erted, respectively, by the “time-symmetric” and “time- 
antisymmetric” pieces of the (regularized) metric pertur¬ 
bation (cf. [23 for a more precise definition). For geodesic 
motion in the equatorial plane of a Kerr black hole, L“ 
can be thought of as a function of only r and r along the 
orbit. The particular time symmetry of such geodesics 
then implies [25] 

^cons(G ^) = \ [^“(l r) + S(„)F“(r, -r)] , (40) 
J^diss(G0 = I [F^(r,r) - s^a)F^r, -r)] (41) 

(no summation over a), where S(i) = —1 = S(^p) and 
S(r) = +1- This gives a simple prescription for construct¬ 
ing and along geodesics, given the full GSF. 

For circular orbits we have F°‘{r,f) = F°‘(r,—f)^ 
meaning L*, are purely dissipative while L” is purely 
conservative. In general, however, each component has 
both dissipative and conservative pieces. Of particular 
interest to us will be nearly-circular orbits with |r| <?; 1. 
Along such orbits we may write, to leading order in |r|, 

F^ans^rF^ir), FS,,,^F^{r) (42) 

for a = t,4>, and 

i"diss=^rF[(r), (43) 


where Eq and L" are some functions of r only. 

Equations (]40])-(43) are applicable, at leading order 


in ? 7 , even for an orbit that is slowly evolving under the 
GSF effect. In that case the GSF depends also on the 
instantaneous self-acceleration, but that dependence ap¬ 
pears only at subleading order in 77 . At leading order, 
Eqs. (40)-(43) maintain their form at each point along 
the orbit. 
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The GSF integrals AE and AL can be related, in cer¬ 
tain situations, to asymptotic fluxes of energy and angu¬ 
lar momentum in gravitational waves. This was estab¬ 
lished rigorously in Ref. m for a trajectory starting and 
ending at infinity]^ A similar balance relation has been 
argued to hold also for adiabatic inspiral orbits around 
a black hole, subject to a suitable averaging over many 
orbital periods [5511221 • In both scenarios, the contribu¬ 
tion from to the integrals AE and AL (taken from 
T = —00 to r = -l-oo) vanishes at leading order, by virtue 
of the orbital symmetry expressed in Eq. (40). This guar¬ 
antees that the radiated fluxes balance the work done by 
the dissipative piece of the self-force alone, as expected. 


C. ADM energy and angular momentum 


Our analysis in the next section will require knowledge 
of the total, conserved ADM energy and angular momen¬ 
tum contents of the spacetime in the above setup. Specif¬ 
ically, we will need expressions for Aadm and Tadm in 
terms of Eoo and L^o (and the two masses, M and /r), 
correct through A subtlety is that ADM quan¬ 

tities are most conveniently evaluated in a “center-of 
mass” system (and, at the required order, would include 
a contribution from the black hole’s “recoil” motion), 
whereas E^o and are components of the particle’s 
four-velocity, defined in a coordinate system centered 
around the black-hole. 

In our setup, the ADM quantities are most easily eval¬ 
uated on a hypersurface of constant t —M, where 
the binary separation is r ^ M. In the limit t —> —oo 
(r —>■ oo), the gravitational interaction energy vanishes 
and does not contribute to Eadm- Working at that 
limit, we assume that, for the purpose of calculating 
ADM quantities, the black hole-particle system may be 
replaced with that of two relativistic pointlike particles 
in flat spacetime. Aadm is then simply the sum of the 
two relativistic energies in the center-of-mass frame, and 
Tadm is similarly the sum of two angular momenta (with 
respect to the center of mass), plus the spin of the black 
hole. 

Appendix [A| gives the details of this calculation, which 
is straightforward. The result is 


Aadm = M 


1 -I- rjE^ 



1 ) 


+ o{v^), 


(44) 


Tadm = M {a + rjL^ - rj'^L^E^) + o{ri^). (45) 


^ The configuration considered in Ref. m had no black hole in 
it, but the authors argue convincingly that a similar conclusion 
would hold also in the black hole case, if fluxes down the event 
horizon were accounted for in the balance equation. 


D. Critical orbits 

In the geodesic case we have introduced the function 
Lc{E), which we now interpret as the critical value of 
the angular momentum for a given energy: Geodesic or¬ 
bits with L > Lc{E) scatter back to infinity, while ones 
with L < Lc{E) fall into the black hole. This type of 
critical behavior carries over to the GSF case, though ra¬ 
diation losses then introduce a subtlety, since orbits that 
are initially scattered may fall into the black hole at a 
subsequent approach. However, we may still speak of a 
critical threshold for an immediate capture, which sepa¬ 
rates (in the space of initial conditions) between orbits 
that scatter at first approach and orbits that do not. A 
detailed analysis of this critical behavior was given in 
Ref. 130] for orbits in Schwarzschild spacetime (working 
in the first-order GSF approximation, as here), and in 
the following discussion we assume the same qualitative 
behavior applies in the Kerr case too. 

In particular, we assume there exists a critical value 
Too = Too,c(Too) that separates between the two pos¬ 
sible outcomes. The initial conditions {Too, Too ,c(Too)} 
thus define a one-parameter family of “critical orbits”. 
Let us denote by Ec(t]Eoo) and Tc(r;Too) the func¬ 
tions E{t) and T(r) corresponding to a critical orbit 
with a given Too [so that Ec{t —>■ —oo;Too) = Too and 
Tc(t —>• —oo;Too) = Too,c(Too)]- Unlike in the geodesic 
case where critical geodesics of different T are disjoint, 
in the GSF case all critical orbits join a global attrac¬ 
tor, which is the perfectly fine-tuned orbit that evolves 
radiatively along the sequence of unstable circular or¬ 
bits starting at the light ring and ending at the ISGO, 
where it plunges into the black hole. Figure 1 in Ref. 
m illustrates the evolution of the critical orbit along 
the attractor, and see also Fig. below. 

Let us dehne the “GSF correction” 



• ^c(j~ ^ ^oo') -^c(^CXd)5 

(46) 

(^Eqq ) 

:= SLc{r -oo;E^). 

(47) 


(5Too is the GSF-induced shift in the critical value of Too 
at a fixed Too. It may also be interpreted in terms of a 
GSF correction to the critical impact parameter. We as¬ 
sume that the difference SLc{t; Too) remains small [0{ri)] 
during the approach, which should be the case in any rea¬ 
sonable gauge. However, clearly, that difference ceases to 
remain small as the critical orbit joins the global attrac¬ 
tor and evolves along it; then the meaning of 6Lc{t', Too) 
as a small GSF correction is lost. 

For our analysis of overspinning orbits in the next sec¬ 
tion, we will require an explicit expression for (5 Too(Too) 
in terms of GSF quantities. It is instructive to derive this 
relation first with the dissipative piece of the GSF turned 
off, i.e. replacing the full GSF with its conservative piece 
(in which case the global attractor disappears, and crit¬ 
ical orbits of different Too remain disjoint). Let us call 
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the resulting quantity 6L'^’^‘^{Eoo)- As a second step we 
will restore dissipation and consider its effect. 


1. Conservative GSF effect 


With dissipation turned off, the critical orbit becomes 
exactly stationary at r —>■ oo, where it joins an un¬ 
stable (nongeodesic) circular orbit of radius R{Eao) = 
R{Eoo) + SR- Here R{Eao) is the geodesic relation given 
in Eq. (14), and SR is a conservative GSF correction. 


To obtain dT™"® we first substitute E and L from Eq. 
( [37| ) into the radial equation of motion ( [3^ , replacing Loo 
with Lc(Loo)+dL“"®(Loo), where Lc{E^ is the geodesic 
relation given in Eq. We then demand dr/dr = 0 as 
well as = 0 at r = .R through 0{ri). At leading 

order in e this yields two algebraic equations for the two 
0{r]) unknowns dL™''® and SR, given E^o and the GSF. 
The solution is 


with a Lyapunov exponent A = M~^{3E‘^ — 1 )^/^ at 
leading order in e (and ignoring the small effect of the 
GSF). The choice Tc = —X~^\ogr], for example, gives 
Sr{Tc) ^ 77, and the quasi-circular piece of does 

not contribute to dL“"® at leading order in r]. 

Our discussion assumes that AE{tc) and AL{tc) 
[hence also AW{tc)\ are 0(77) quantities, i.e. that the ac¬ 
cumulated GSF-induced positional shift in the orbit dur¬ 
ing the approach is a small, 0(77) quantity. This should 
be the case in any reasonable gauge. Under this assump¬ 
tion, we note that the value of the integral AW remains 
unchanged, at leading order in 77, if in Eq. ( 50 ) we replace 
the integration along the actual, GSF-perturbed orbit, 
with an integration along the critical geodesic of energy 
Loo. This can be exploited to simplify actual calcula¬ 
tions: To compute dL“"® at leading order in 77 requires 
only an evaluation of the GSF along a fixed geodesic, and 
there is no need to consider the back-reaction from the 
GSF on the orbital trajectory. 


(5L“""(Loo) = 2MA£;“""(oo) - AL'=™"(oo), (48) 


and SR{E^) = 0{e)0{r]). Here AL™''" and AL“"" 
are the same as AE and AL of Eq. ( [M| , but with 
Fa —>■ L™"®, and with the GSF integrals evaluated along 
the critical orbit with energy-at-infinity E^o ■ The precise 
dependence of SR on the GSF will not be needed, but we 
note that the 0{erj) GSF correction to the radius of the 
critical circular orbit is reassuringly small compared to 
the 0 (e) radial distance to the light ring. 

To simplify the appearance of subsequent equations, 
let us from now on use units in which M = 1. This, in 
particular, makes our “tilde” notation redundant (with 
L = L, etc.) and 77 becomes interchangeable with 77 . 
Recalling our W notation from Eq. (21), we rewrite Eq. 
(IdSl) as 


dL™"^(L;oo) = AW™"'^(oo), 
where AVF“°"® represents the conservative piece of 


(49) 


AVF(t) := 2AE{t) - AL(t) 

= - 7 ?-^/’ {2Et + E^)dT. ( 50 ) 

J —00 

The quantity AVF“"®(oo) is the limit t —)■ 00 of 
AVF“"®(r). Does this limit actually exist? The answer is 
positive, since both E™"® and L™''® vanish exponentially 
fast in T as the orbit approaches the limiting circular or¬ 
bit at r —>■ 00. 

To make this last statement more precise, let us split 
the r integral into an “approach” piece, and a 

“quasi-circular” piece, with Tc chosen so that Sr(Tc), 

where Sr(T) := rij) — R, is already very small. For a 
small Sr we have the form ~ rEit(r) [Eq. (42)] and 


similarly for Thus E^°'^^dT ~ —Fit{R)Sr{Tc), 

and sim ilarly for the (j) component. A local analysis of 


Eq. (35) near the limiting circular orbit gives Sr e 


— At 


2 . Full GSF effect 


Now restore dissipation. The fine-tuned critical orbit 
no longer settles into a strictly stationary motion, but 
rather it continues to evolve radiatively, in an adiabatic 
fashion, through a sequence of unstable circular orbits 
of decreasing energies (hence increasing radii). With a 
perfect fine-tuning, the orbit can reach the ISGO be¬ 
fore plunging into the black hole—a scenario illustrated 
in Fig. [^ A relation between the degree of fine-tuning 
and the amount of energy loss was derived in Ref. EH] 
(for the Schwarzschild case): Rewriting their Eq. (124) 
in terms of angular momentum, we have the scaling 
5L/Loo,c exp[(Lf — Li)/ 77 ], where SL (not to be con¬ 
fused with the GSF shift SLao) is any small perturbation 
in the value of the initial angular momentum off the crit¬ 
ical value Loo c) and E{ — L, is the resulting change in 
specific energy as the orbit progresses along the global 
attractor. To achieve an 0(1) change in the specific 
energy requires an “exponentially delicate” fine-tuning, 
5L/Loo,c ~ exp(-l/77). 

For our analysis we do not require knowledge of the 
perfectly fine-tuned angular momentum at that level. We 
need Loo,c through 0 ( 77 ) only. Fine-tuning at 0 ( 77 ) [cor¬ 
responding to SL = 0(77)] guarantees only Ef — Ei = 
0{r] In 77). Therefore, for the purpose of determining Lao,c 
through 0 ( 77 ), it is sufficient to restrict attention to the 
early part of the critical orbit, where AE and AL (spe¬ 
cific values) are still 0(77 In 77) at most, and have not yet 
accumulated 0(1) changes. This observation assists us 
in Appendix where we derive an expression for the 
leading-order full-GSF correction SLadEao)- 

Our main result in Appendix [^ is that 


SLaa{Eaa) = AW{tc) + 0{e)0{r]lnri), (51) 

in analogy with the “no dissipation” case, Eq. ( |4^. Here, 
AW{tc) is the full-GSF integral shown in Eq. (50), eval- 
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FIG. 4. Schematic illustration of the evolution of orbital 
energy along a perfectly fine-tuned critical orbit (solid blue 
curve). The orbit approaches from infinity, becomes trapped 
on an unstable circular orbit, and then evolves adiabatically 
in a quasi-circular fashion before transiting to plunge around 
the ISCO location. Radiative losses are small during the ap¬ 
proach and plunge, but, through fine-tuning, the orbit can 
be made to lose “all” its energy during the quasi-circular 
stage. Intermediate values of the final energy Ef may also 
be obtained by fine-tuning (dashed magenta line). Note the 
orbital radius increases through radiation losses during the 
quasi-circular stage. In the near-extremal case, e <C 1, the 
quasi-circular evolution and final plunge occur within a small 
range of coordinate radii, Ar — 0{e). 


uated along the orbit from infinity and up to the “end 
of approach” time Tc, when the orbit settles into a quasi¬ 
circular motion. Crucially, the contribution to AW from 
the quasi-circular part of the orbit is suppressed by a 
factor of e, so that the precise choice of Tc does not af¬ 
fect the value of AW(tc) at leading order. This assumes 
only that e <C | In 7^1“^, so that the error terms in Eq. 
(51) are negligible compared to AW{tc) = 0(e°)0(?7). 
AU we require of Tc is to be sufficiently late that |r| is 
already very small [specifically, r(Tc) = 0{r])], but suffi¬ 
ciently early that AE is 0(7? In 77) at most. In practice, 
AW may again be evaluated along the critical geodesic 
of energy-at-infinity with the integral in Eq. (51) 


truncated after, say, 4-5 orbital revolutions. Truncat¬ 
ing instead after (e.g.) 10 revolutions should change the 
result by a negligible amount of only 0{e)0{r]). 

The dissipative piece of AW{tc) [call it de¬ 
fined by replacing Fa —f in Eq. (51)] may be ex¬ 

pressed in terms of radiated quantities. Let i5(apr) and 
'C(apr) be the total energy and angular momentum in 
gravitational waves radiated out to infinity and down the 
black hole during the approach. We shall assume that the 


balance relatior FFl 


^^(apr) • ^^(apr) -^(apr) 


= -7?AlE‘^‘""(rc) (52) 


holds at leading order in 77 and in e. Equations ( |49| ) and 
(51) then lead to 


5Laa = - W(,pr)/77, 


(53) 


where subleading terms have been omitted. This reex¬ 
presses 5Lao as a sum of conservative and radiative con¬ 
tributions, the motivation for which will become clear in 
the next section. 

Finally, let us further write £(apr) = '^^Jpr) + ^(apr) 
similarly for £(apr) and ^(apr)! where hereafter super¬ 
scripts ‘-b’ and ’ denote contributions from fluxes to 
infinity and down the black hole, respectively. The fol¬ 
lowing argument, based on the first law of black hole 
mechanics, suggests that must vanish in the limit 

e —>■ 0. If we assume the black hole is not overspun dur¬ 
ing the approach, its horizon’s area should increase by an 
amount 5A satisfying 

= '^(apr) - ^^^'^(apr)- (^4) 


where n = e/-\/2 -|- O(e^) is the horizon’s surface gravity, 
and fin its angular velocity. Since fin = ^ + 0(e), we 
identify the right-hand side of (54) as at lead¬ 

ing order in e. We thus have, m the extremal limit, 
~ e(ci(5A-|-C2A[)^p,.)), with ci,C 2 certain numerical 
coefficients. Since 5A and must remain bounded 

even in the extremal limit, we conclude that Wr , van- 
ishes in that limit. Thus, at leading order in e, Eq. (53) 
becomes 


-W+p^)/77, (55) 

which now features only outgoing fluxes. 

With this we have completed the necessary ground¬ 
work for our overspinning analysis, to be presented next. 


^ The balance | |52[ | does not follow directly from the theorem of 
Ref. [27) . because the approach part of the critical orbit does 
not end at infinity. It may be possible to construct a proof by 
considering a small outward deformation of the orbit (such that 
the new orbit starts at infinity and scatters back to infinity), 
then invoking the approximate symmetry about the periapsis, 
together with the e-suppression of the quasi-circular contribution 
to We shall not endeavour to provide the details of such 

a proof here. 

^ One cannot expect to be able to similarly balance £’(apr) 
£(apr) individually, because the dissipative pieces of HEirc) and 
AL(tc), unlike AW^i""(rc), are sensitive to the choice of Tc al¬ 
ready at leading order. However, such individual balance rela¬ 
tions will not be needed in our analysis. 
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IV. OVERSPINNING WITH THE FULL 
SELF-FORCE 


A. General form of the censorship condition and 
reduction to near-critical orbits 


Starting with a near-extremal Kerr geometry with 
a/M = 1 — e^, consider a particle sent in from infinity 
with specific energy Eco and specific angular-momentum 
Loo at t —>■ —00. The ADM mass and angular momentum 
of the spacetime are given in Eqs. ( 44 ) and ( 45 ) through 
O(n^). We assume the particle crosses the event hori- 
zorjj at some retarded time Uh, and we let and £■*" 
be the total energy and angular momentum radiated to 
null-infinity up until Uh (with u/j —> oo if the post-capture 
geometry relaxes to a Kerr black hole). Then the Bondi 
mass and angular-momentum of spacetime at retarded 
time Uh are Ladm — and Ladm — respectively. 
Overspinning is avoided if and only iQ 


(^^ADM — Y ^ ^ADM — ■ 


( 56 ) 


To rule out the overspinning scenario, we need to show 
that this inequality holds for all Eao,Lao and for all rj, e. 
Since we work in first-order perturbation theory, we only 
demand that ( 56 ) is satisfied at leading order in 77. We 


also assume e <C 1 and keep only leading terms in e, but 
we do not a priori restrict the magnitude of e relative to 
that of r]. We shall refer to the inequality ( 56 ) as the 
censorship condition^ 


Substituting from Eqs. ( | 44 | ) and ( [ 4 ^ , the censor¬ 
ship condition becomes 


+ pWoo + ??^(1 + LooE^o — L^) 

+ (77^00-L+)"-W+> 0 , ( 57 ) 

where Woo '■= 2£’oo —Loo, W"'" := 2L+ —£+, and we have 
omitted subleading terms of 0(77^). Note how the various 
terms here scale with 77. The quantities E^o and Lqo 
( hence also Woo) are specific values, thus a priori they 
are 0 {rf). The radiated energy L+ is generically 
but may accumulate at 0(77) for orbits that are fine-tuned 
to evolve along the global attractor; it is to allow for 


such orbits that we have kept the terms 2 rjEoo£^ and 
{£'^Y in Eq. ( (^ . The quantity is likewise 0 {rf) 
generically and up to 0 (rf) with fine-tuning, but, as will 
be shown below, in the latter case the 0(77) term is also 
proportional to e. 


Inspecting Eq. ( 57 ), we observe that, for all captured 


orbits that are not sufficiently close to criticality, the term 
jyWoo is 0(77) and positive, so the censorship condition 
( 57 ) is trivially satisfied at leading order in 77 and e. Vio¬ 


lation of ( 57 ) (hence overspinning) may only be achieved, 
potentially, if L^o is tuned so that = 2 Eao + 0(e, 77), 
giving Woo = 0(6,77). It is therefore sufficient to restrict 
attention to this class of orbits, to be referred to in what 
follows as “near-critical”. Formulated on near-critical or¬ 
bits, the censorship condition takes the sufficient form 


e" + 77W00 + 77^(1 + L^) + ipE^ -£+y-W+> 0 , ( 58 ) 


where we have dropped 0 {erf', rf) terms. This is required 
to hold for each member of the reduced two-parameter 
family {Loo, Loo} with - 2 E^ = 0 {e,r]). 

To proceed, we need to make more precise the dis¬ 
tinction between “fine-tuned” orbits and “generic” near- 
critical orbits that are not fine-tuned. Referring to 
Fig. i, let Ef he the final value of the specific energy 
with which the particle plunges into the hole; and let 
Loo,c(Loo) be the perfectly fine-tuned value of Loo, for 
which the orbit joins the global attractor and evolves 
along it to the ISCO. Assuming the universal scaling 
Loo Loo,C (Loo) '^ooi)/ 77 ] [ 30 j . near critical 

orbits as defined above generically have Ef — Loo = 
0(77ln?7e) [here we neglect the 0(77) difference between 
Loo and L^]. Calibrating Lao at higher order in rj,e [so 
that Lao —Lao,c = 0(77”, 6^) with some n,k > 1 ] does not 
qualitatively change this generic scaling of L/ — Loo. To 
achieve Ef — Loo = 0(1) requires an exponentially accu¬ 
rate tuning, i.e. Loo—Lao,c ~ exp(—1/77). In what follows 
we use the 77 scaling of L/ — Loo to distinguish between 
generic and fine-tuned members of the near-critical fam¬ 
ily: The former admit Ef — Loo = 0(77 In 775), and the 
latter Ef — Loo = 0 ( 1 ). This distinction can also be 
formulated in terms of the radiated quantities £^ or 

0(77^ In776) (“generic”), ( 59 ) 

0(77) (“fine-tuned”). ( 60 ) 


® More pedantically, we refer here to the crossing of a marginally 
trapped surface; spacetime need not contain a global horizon. 

^ We do not know, and for our purpose do not need to know, the 
future evolution of spacetime beyond retarded time in the 
hypothetical case where ( |56| l is not satisfied. The likely scenario 
involves the formation of a naked singularity and a breakdown 
of predictability for u > Uh If ( |56[ | is satisfied, then, by 
“no-hair” theorems, ge ome try should relax to a Kerr black hole. 

® It may be argued that ( |56| l is guaranteed to hold (with a strong 
inequality) by virtue of the third law of black-hole mechanics 
m, though it is clear that some of the third-law’s assumptions 
are not satisfied within our model—see [9] for a discussion. Even 
if that can be established, it is still of interest to explore the 
physical mechanism that enforces the third law in our setup, 
which is what our study aims to achieve. 


B. Further reduction to critical orbits 


The inequality ( 58 ) is still a condition on a two- 


parameter family of orbits. Ignoring fine-tuned orbits for 
now, it is possible—and beneficial—to reduce it further 
to a sufficient condition formulated on a one-parameter 
family. We achieve this by minimizing the left-hand side 
of Eq. ( 58 ) over all near-critical orbits for a given Eao ■ 


We argue that the minimizing orbit is one with Lao tuned 
to its critical value at least through 0(77,6), namely 


Loo = 2 Loo + 6 a/ 6 L^^^-|- 5 Loo(Loo) + o(?7,6), ( 61 ) 
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where we have recalled Eq. (16), and 5Eoo(Eoo) is the 


0{ri) GSF term derived in the previous section. To see 
this, note first that Woo = — Loo is trivially mini¬ 

mized by Loo = Loo,c{Eoo), since Too,c(E'oo) maximizes 
Loo (over all captured orbits of a fixed Eoo) by definition 
of a critical orbit. This means that, to minimize Woo 
through 0{r],e) (higher orders are irrelevant in our ap¬ 
proximation) it suffices to demand Loo — Too,c(Eoo) = 
o{r],e). Then also note that the two radiative terms 


on the left-hand side of (58) are insensitive, at relevant 
order, to variations of Too within the family of nearly- 
critical orbits for a fixed Eoo- For generic orbits, the 
term {rjEoo — is simply rfE^o at leading order, re¬ 


calling Eq. (59). As for the term —W+, we note that 


the contribution to that term from the approach part of 
the orbit, which is already O(ry^), is not sensitive, at that 
order, to 0{r],e) variations in Loo- Meanwhile, the con¬ 
tribution to yV~^ from the quasi-circular part of the orbit 
is of 0{e)0{r]‘^ In rje) at most (the occurrence of the fac¬ 
tor e will be explained below) and hence negligible in Eq. 
(58), assuming only e |ln 7 y|“^. 

Thus, discounting fine-tuned orbits, we find that the 
entire left-hand side of Eq. (|58[) is minimized by Loo as 


given in Eq. (61). A new sufficient version of the censor¬ 


ship condition may therefore be written as 

- T^ey'QEl - 2 -r,SLoo+ + E^) 

+ (??Eoo-T+)"-W+> 0 , (62) 

which, at the relevant, leading order, is a condition 
on the one-parameter family of (generic) critical orbits 
parametrized by Eoo alone. 


It should now be noted that the condition (62) also ap¬ 


plies to fine-tuned orbits [whether or not they minimize 
the left-hand size of (58)], simply because such orbits al¬ 
ways satisfy Eq. (61). However, for fine-tuned orbits the 


condition still involves two parameters, conveniently cho¬ 
sen as Eoo and Ef. Different values of E/ correspond to 
a fine-tuning of Loo at an exponential level. In princi¬ 
ple, any value of Ef in the range Eisco ^ Ef ^ Eoo may 
be obtained this way. To rule out overspinning by fine- 
tuned orbits, the censorship condition (621 must hold for 


all {Eoo,Ef} with Ef in the above range. 
Observe that in Eq. (62) we have 


= vl: 


cons 

oo 




-t 

(apr) 


;+ 

(end)’ 

(63) 

f g’ 



where in the first line we have recalled Eq. (|M), VV(q(, 

, ong the 

quasi-circular part of the orbit, and is the piece 

coming from the transition to a final plunge into the black 
hole and from the plunge itself. It follows that only the 
conservative piece of the shift 5Loo actually enters our 
condition: 

- r^e^&Elo - 2 - ??5T“- + rf[l + Elo) 

+ (^Eo.-T+)^-W+,)-W+„,)>0. (64) 


In this last form, conservative and dissipative terms 
of the GSF feature separately. The former are associated 
with the approach leg of the orbit, and the latter accumu¬ 
late during the adiabatic evolution along the attractor. 
In Appendix [O we combine results by Ori and Thorne 
[32], Kesden |33| and Mino and Brink |34|, to argue that 
the term is always subdominant and negligible in 

Eq. (64). We shall therefore omit that term in the rest 


of our discussion. 

In subsection IIVDI below we will show that the ra¬ 
diative term scales as 0{e)0[r]{Eoo — Ef)]. This 

term can thus feature at leading order in Eq. (64) only 
for fine-tuned orbits, for which Eoo — Ef = 0(ly Like¬ 
wise, terms involving T"*" feature only for fine-tuned or¬ 
bits and are negligible otherwise. On the other hand, 
the conservative term 77 T™''® is always featuring 

in the censorship condition regardless of fine-tuning. An 
important consequence is that dissipative effects of the 
GSF enter the censorship condition only for fine-tuned 
orbits. This seems consistent with suggestions made in 
earlier analyses [131 ESI |35| (in which fine-tuning has not 
been considered). 


Below we further simplify the condition (64), and re¬ 
formulate it explicitly in terms of Eoo alone (for generic 
orbits) or Eoo and Ef alone (for fine-tuned ones), with¬ 
out reference to p and e. We shall consider separately 
the cases of generic and fine-tuned orbits, starting with 
the former, simpler case. 


C. Censorship condition for generic orbits 

As mentioned above (and shown in the next subsec¬ 
tion), without fine-tuning the radi ativ e terms and 


^(qc) become subdominant in Eq. (|64|) and drop out of 


it. The censorship condition then reduces to 


with 


-I- peF p^H > 0, 


E:=-V6E^-2, 


H-.= l + 2Eio- 5LZ^ 


Here we have made the ry-scaling of 5L\ 
introducing the shift-per-eta 


cons 

00 


5T“- :=ry-ijT™“, 


(65) 

( 66 ) 

(67) 
explicit by 

( 68 ) 


which should have a finite (nonzero) limit 77 —> 0. 

For the overspinning scenario to be ruled out, the in¬ 
equality (65) must be satisfied for all p,€ > 0 and all 
Eoo > 1- The condition can be written in the equivalent 
form $ := -I- aE H > 0, with a := e/p > 0. At 

fixed Eoo, 4> is quadratic in a, with a minimum value 
$(q; = —E/2) = F[ — E"^/L To guarantee $ > 0 for all 
Eoo and all a > 0 (hence all 77 , e > 0 ) we must demand 
F[ > E'^/A] if < E^/4 for some Eoo, then for that Eoo 
there exist 77 , e values for which overspinning is achieved. 
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In that way, H > F"^/A is both sufficient and necessary 
for overspinning to be avoided. Inserting the values of F 
and F[, the censorship condition takes the simple form 


SL": 




3). 


(69) 


Overspinning is averted (for orbits that are not fine- 
tuned) if and only if (691 is satisfied for each member of 


the one-parameter family of critical orbits with > 1 , 
in the limit r],e ^ 0. 

Equation (69) states our final result for generic orbits. 


As already mentioned, it involves only conservative GSF 
effects, specifically the shift in the critical value of the 
angular-momentum-at-infinity (at fixed Eao) due to the 
conservative piece of the GSF. For easy reference, we 
give here the explicit formula for in terms of GSF 

components: 


5L—(Fo.) = hm 4 

£->■0 J 


pOO 

Jr, 


(70) 

where we have recalled Eqs. (49) and ([50|. The integra¬ 


tion is carried out along the critical geodesic of specific 
energy Eoo on a background with spin a/M = 1 — e^, 
starting at the unstable circular orbit of radius = 
R,:{e,Eao) and ending at infinity. 

Inspecting Eq. (69), it may seem peculiar that over¬ 


spinning may be averted even for some positive values of 
A positive would seem to mean that the 

GSF increases the critical impact parameter, allowing in 
particles that would otherwise be scattered away. How¬ 
ever, we must recall that the shift is defined not 

with respect to the physical, ADM angular momentum, 
but with respect to the quantity L = which (while 
convenient to work with in practice) does not have a clear 
invariant meaning beyond the geodesic approximation. 
To rewrite (69) in a more physically insightful way, let 


us, then, recast it in terms of ADM quantities, as follows. 

First, let us introduce the specific quantities and 

L^dm defined through 


^-^ADM •— A^ADM Af, 
M^adm •“ Aadm Ma, 


(71) 


which may be thought of as the particle’s contributions 
to the total ADM energy and angular momentum of the 
system. Then, denote by <5 Tadm('®adm) ^de shift, due 
to the conservative GSF, in the critical value of Tadm. 


a fixed ^ short calculation, based on Eqs. (|44|) 


and ( |45[ ), gives 

SLluuiEluM) = SL‘^°^^iE^)-v{El + l)+0{v^). (72) 

Thus, in terms of <5 Tadm '^^adM’ ^de censorship 

condition (69) becomes 


'^^ADM — 2 


^ ) 
00/5 


(73) 


where on the right-hand side E^ may be replaced with 
^ADM relevant order. 


The alternative form (73) is now more intuitive: For 
unbound orbits {Eoo > IjTthe GSF averts overspinning 
if it shifts the critical value of the (ADM-related) angu¬ 
lar momentum by a sufficiently negative amount, which 
depends only on E^o- In the marginal case of E^o = 1 
(where overspinning is marginally prevented already in 
the geodesic case), the shift i^Tadm need only be non¬ 
positive. We are not aware of any a priori argument to 
suggest that (^Tadm niust necessarily be nonpositive for 
all Eoo > 1. Verifying this would need to await a nu¬ 
merical calculation. Any counterexample would imply a 
direct violation of cosmic censorship. 

Let us make a few more points about the condition 
(69). First, in its form (73) it is manifestly gauge in¬ 


variant (within a class of suitable asymptotically flat 
gauges) despite the gauge dependence of the local GSF 
featuring in [Eq. (70)]. The condition involves 

only quantities that are defined and evaluated at infinite 
separation, namely the specific energy E^o (or T^adm) 
and angular-momentum shift each having a clear 

gauge-invariant physical meaning. The evident invari¬ 
ance of our final condition is reassuring. 


Second, as already mentioned, the condition that (69) 
is satisfied for all Eoo > 1 is both sufficient and necessary 
for overspinning to be avoided within the scenario con¬ 
sidered here. It is a sujficient condition only in the sense 
that it guarantees no overspinning occurs for sufficiently 
small mass-ratio 77 ; since we work in the first-order self¬ 
force approximation, we cannot make the statement any 
stronger. Equation (69) describes a necessary condition 
in the sense that its violation for any E^o would mean 
there exist (small) 77 values for which overspinning oc¬ 
curs. 


Finally, the condition (69) involves the single param¬ 
eter Eoo, and the task of testing whether it is satisfied 
amounts to evaluating a single function of E^o, namely 
5L^ff^{Eao)- The perturbative parameters themselves, 77 
and e, do not feature in the final condition. This is ex¬ 
pected, given our first-order perturbative treatment and 
the fact that GSF effects (including ADM terms) appear 
in the overspinning condition already at leading order. 
It is precisely because of this “order mixing” that one 
cannot neglect the GSF in considering the overspinning 
problem, and why there is no sense in which the geodesic 
limit may be said to provide a useful approximation here. 


D. Censorship condition for fine-tuned orbits 

It is not a priori clear whether fine-tuning favours 
the overspinning scenario or disfavours it: The answer 
depends on the details of the radiative evolution along 
the attractor. Indeed, for fine-tuned orbits the radiative 
terms and feature already at leading order in 

Eq. (64), and cannot be neglected. We may again write 
Eq. (64) in the form (65), with F and FI replaced with. 
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respectively, 

H = l + El^ + {E^-£+)^ (74) 
Here we have introduced the rescaled quantities 

:= := (75) 

which should have finite (nonzero) limits e, r; —>■ 0 for fine- 
tuned orbits [that = 0{e)0{r]) will be discussed in 

the next two paragraphs]. 

It will prove beneficial to reexpress E and iJ in terms of 
the absorption-related quantities £~ and place 

of £+ and W(q(.)- This is easily done for H, noting E^o — 

£'^ = Ef + £~ at the relevant, leading order. As for F, 
we start by writing = W(qc) — , where, under 

the assumption of adiabaticity, the total W(qc) may be 
expressed as an integral over the local GSF: 


W(qc) = 


(2^diss p. 


idiss^ 


dE/E. 


(76) 


Here we have used Eq. (50), changing the integration 


variable from r to specific energy E, and assumed a bal¬ 
ance relation as in Eq. (52). We have also neglected the 
subdominant [0{r]‘^)] amount of radiated energy during 
the approach, replacing the initial energy of the quasi¬ 
circular motion with E^ . Then, following the method of 
Appendix[cf. Eq. (^)], we use = 0 to obtain 


^disS-^^diss/2 ^ _(yr/yt)^diss_ 


V6E2 _ 2 ‘ 


diss 


, (77) 


where subdominant terms in e have been omitted. The 
contribution from the term oc to the integral in (76) 
can be evaluated following the same steps as in Appendix 
[b] [see the paragraph containing Eq. (B5)], and shown 
to be of only 0{e)0{rf') (or smaller) —hence negligible. 
The contribution from the term oc can be evaluated 
explicitly upon replacing Ft = fJ-E, giving 


and Ef only. While £'^ is necessarily positive, the ab¬ 
sorbed energy £~ may be either positive, or—due to 
superradiance—negative, depending on E^o and Ef. Cir¬ 
cular equatorial geodesics are superradiant for 17 < fin, 
which, in the extremal limit, translates to F < Thus, 

is necessarily negative for any Ef < E^o < The 

sign (and magnitude) of £~ for other values of Eao and 
Ef, as well as the sign (and magnitude) of remain 

to be determined numerically. The quantity on 

the other hand, is easily shown to be negative definite. 
In fact, Eqs. ([84|) and (| 8 ^, given below, imply 




(80) 


Note this means that F in Eq. (74) may change sign, 
depending on E^o, Ef. 


To proceed, we once again write the condition (65) (for 
the barred quantities) in the form $ := + aF + H > 0, 
with a = t/rj. Here, however, the sign of F is not known a 
priori, which somewhat complicates matters. For F < 0, 
$ has its minimum at i>(a = —F/2) = H — j^, so 

the condition becomes H > F^/4 as before. However, 
for F > 0 the condition $ > 0 is satisfied trivially for all 
H > 0, and violated trivially for all F < 0 (by choosing a 
sufficiently small a). In that case, therefore, a necessary 
and sufficient condition for $ > 0 to hold for any 77 , e is 
F > 0. In summary, we obtain 


F > (min{F/2, 0})^ 


(81) 


as a necessary and sufficient condition for overspinning 
to be averted for all 77 , e. In this condition, F and F are 
both functions of the two independent parameters E^o 
and Ef. To rule out overspinning we must require that 
(81) is satisfied for all E^ > Ef{> Fisco)- 
^Evaluation of the condition ( [M| ) requires knowledge 
of the radiative quantities £^ and (^’^ addition to 

To conclude our discussion, we now give conve¬ 
nient expressions for these two quantities in terms of a 
single function of one variable, namely the ratio 


W(qe) = -776 (V6F^ - 2 - ^6F2 - 2) . (78) 


Thus, Eqs. (74) are obtained in their alternative form 


F = -,/ 6 F| - 2- 


W, 


(qc) 

-^2 


H = l + E^ + {Ef+£-y -SE 


(79) 


We note that Eq. (78) establishes the scaling kV(qc) = 
0{e)0{ri) for fine-tunea orbits. The first-law argument 
used in the pre vious section [refer to the discussion 
around Eq. (54)] can also be used to show = 

0{e)0{r]). This then establishes the scaling Wjq;.) = 
0{e)0{f]) assumed above. 


In both Eqs. (74) and (79), the radiative quantities 
£^ and should be thought of as functions of E^o 


n[E) := 


f+(F)’ 


(82) 


where £'^/~ (F) are the outgoing/incoming fluxes of en¬ 
ergy in gravitational waves sourced by a particle on a 
circular geodesic, evaluated in the extremal limit at a 
fixed specific energy F. [In deviation from our notation 
elsewhere, here and in the next two paragraphs an over¬ 
dot denotes differentiation with respect to (any suitable) 
coordinate time.] We note 77. < 0 for F < the super¬ 
radiance regime in the extremal limit. 

First, we use the specific energy F as a parameter along 
the global attractor, to write 


„ r^i p+ 

£+ = / -^dE = - 
77 F 


£+ 

Ie,^ £+ + £- 


dE, 


(83) 
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where we assumed the balance relation rjE = 
applies during the adiabatic evolution along the attrac¬ 
tor. Thus, 


£+{E^,Ef) 


dE 

i + n{E)^ 


and, similarly, 


£ {Eao,Ef) 



nE) 

iEn{E) 


dE, 


(84) 


(85) 


which should be evaluated in the extremal limit, e —>■ 0. 
Note f ^ 0 in the extremal limit m, while the ratio TZ 

admits a finite, nonzero limit [361137] . Thus, by writing 


£^ as in Eqs. ( |84[ ) and ( |85[ ) we have made it possible for 
the limit e —)■ 0 to be taken before the integration, which 
is advantageous in practice. 

As for we start by writing 


:= ± = -eb{E)£^ 


( 86 ) 


where are the angular-momentum fluxes correspond¬ 
ing to and 


b{E) := 


6E 


y/6E^ - 2 ■ 


(87) 


To derive the second equality in (86), which is valid to 


leading order in e, we have used the small-e expansion of 
the orbital angular velocity at fixed E, 


n=l-\biE)e + 0{e% 


( 88 ) 


together with the general relation £^ = applica¬ 

ble to the radiation from any circular orbit [38] . Thus, 
proceeding as with we obtain 



r^f w+ 

lim / --^dE = 

^rjE 



1 + n{E) 


dE, 


(89) 


and, similarly, 


w-,= EmmiE. 


(qc) 


l + 7^(^;) 


(90) 


Equations (84), (85), (89) and (90) express £^ and 
in terms oiThe single iunction7c(E), left to be deter¬ 
mined numerically. 

In Sec. IVII we will assess the numerical task of eval¬ 
uating our censorship conditions, review the status of 
relevant existing GSF codes, and comment on how they 
would need to be modified in order to provide the neces¬ 
sary data. But first, in the next section, we take a short 
detour to explore an alternative approach to the deter¬ 
mination of which offers a practical advantage. 


V. REFORMULATION IN TERMS OF 
REDSHIFT VARIABLE 


Our final overspinning conditions (69) and (81) feature 
the critical angular-momentum shift whose eval¬ 

uation, through equation (70), requires an integration of 
the GSF from infinity along critical geodesics. As we dis¬ 
cuss in the next section, this step is the main stumbling 
block when it comes to evaluating the conditions using 
currently available GSF codes. The integration from in¬ 
finity comes about, essentially, because of the need to 
relate the local properties E and L of the particle just 
before it falls into the black hole, to ADM properties of 
spacetime defined at infinity. This would have been un¬ 
necessary if we had available explicit formulas for Fadm 
and Ladm (or for the corresponding Bondi quantities 
7?adm — £^ and Ladm — >C+), correct through 0(E), 
for the configuration of a particle in an unstable circular 
orbit around a Kerr black hole. Furthermore, given such 
formulas we would have been able to relax the require¬ 
ment that the particle is sent in from infinity, and explore 
the possibility of overspinning with “bound” orbits. (We 
recall our result from Sec.|n]that bound geodesics cannot 
overspin; however, in principle, there remains the possi¬ 
bility that GSF effects change this situation.) 


By good fortune, suitable formulas have been proposed 
very recently, in Ref. m- The expressions, to be pre¬ 
sented below, were obtained using (and in agreement be¬ 
tween) two independent frameworks. One is the Hamil¬ 
tonian approach of Isoyama and collaborators [39] . in 
which the conservative portion of the orbital dynamics 
is described (through first order in rj beyond the geodesic 
approximation) in terms of geodesic motion in a cer¬ 
tain effective smooth spacetime. The other is based on 
the recently proposed “first law of binary black-hole me¬ 
chanics” [40l [41] (itself a limiting case of the generalized 
law established in [42|), which relates ADM properties of 
a helically-symmetric binary system of post-Newtonian 
particles to the so-called “redshift” of the particles (see 
below). Neither frameworks is a priori guaranteed to 
correctly describe the strong-field dynamics in the black- 
hole-particle system relevant to us, but some evidence 
suggests that they might (we return to discuss this point 
at the end of the section). 


The said results, as they are stated in m, apply to 
a particle in a circular equatorial orbit (stable or un¬ 
stable) around a Kerr black hole, ignoring the dissipa¬ 
tive piece of the gravitational interaction (or, more pre¬ 
cisely, time-symmetrizing the gravitational perturbation, 
so that spacetime admits a global helical symmetry). 
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They express the Bondj^ energy and angular momen¬ 
tum of that configuration, through 0 {rf')^ in terms of 
Detweiler’s redshift variable [ 55 ] 

^ := {uT\ ( 91 ) 

where u* is the t component of the four-velocity on the 
circular orbit, and overhats, recall, denote properties of 
the GSF-corrected orbit. The usefulness of such relations 
is in the fact that a computation of z requires only GSF 
information for circular orbits, and there is no need to 
integrate from infinity. Such information is essentially 
accessible to existing GSF codes. 

Following [TSj , let us formally expand the redshift z in 
powers of rj, in the form 

z = zoi^)+r]Zi{D,) + 0{r]'^), (92) 


where ri(= d(j)/dt) is the circular orbit’s angular velocity. 


zo = (1 - aO)i/2 [i + _ 3(ikm)2/3(i _ af])i/3 


1/2 


(93) 

is the geodesic limit of z, and r]Zi{n) is the 0{f]) GSF 
correction, defined for a fixed value of fl. According to 
Ref. |16j . the Bondi energy and angular momentum of 
the circular-orbit binary are given, through 0 ( 77 ^), by 


= M + fiEl , = Ma + , (94) 


where 






dz 


with 


z(n) = zo(fl) -H ^rizi{n) + 0{if). 


(95) 


(96) 


The label ‘sym’ is to remind us that these Bondi prop¬ 
erties are defined in a time-symmetrized (“half-retarded- 
plus-half-advanced”) spacetime. The function Zi(f2) ex¬ 
plicitly determines and Tg™ through 0{ri^). 

We are now reaching the crux of our discussion. Con¬ 
sider a critical orbit, subject to the conservative GSF 
alone (dissipation ignored), which asymptotes to a cer¬ 
tain unstable circular orbit at t —> 00 . Let E^™{u) 
and L^™(m) be the Bondi energy and angular momen¬ 
tum of the corresponding time-symmetrized spacetime, 
with u a suitable retarded-time coordinate. At u —>■ 00 , 
these quantities must approach the corresponding Bondi 


® First-law literature |40l 1411148| usually alludes to ADM proper¬ 
ties, which are defined even in helical symmetry within the PN 
context in which these works operate. In the context of black hole 
perturbation theory, the first-law results should be interpreted as 
referring to Bondi properties. See also mi. where first-law re¬ 
sults are formulated directly in terms of Bondi quantities for a 
black-hole-particle system. 


quantities of the asymptotic circular-orbit configuration, 
as given in Eq. (94). Furthermore, 


E^^'^u oo) = Aadm, A§'“(u -)■ oo) = Ladm, 

(97) 

where on the right-hand side we have the ADM prop¬ 
erties of the physical (“retarded”) critical-orbit space- 
time. [That this must be the case follows from E^’^{u —>■ 
oo) = A'ai^m - = -E^adm, where 

and E~ are the total energies flowing, respectively, 
outward at future null-infinity and inward at past null- 
infinity, in the time-symmetrized setup where = E~. 
A similar argument applies to the angular momentum.] 
As a result, we can write Aadm = M + M-^adm 
Tadm = M + /tLadm N in Eqs. ([^], with 


mP _ z n 


rP 

-^ADM 


dz 

dfl 


(98) 


These expressions relate the ADM properties of the phys¬ 
ical critical-orbit configuration to the redshift of the 
asymptotic circular orbit when dissipation is ignored. 

The conservative GSF shift < 5 L^j 3 ]y[(i?^j 3 ]y[) [recall Eq. 
( |72| )] may now be obtained simply by considering the 
0{^) piece of L^Qj^_in Eq. (|^, for a fixed 


pP 

-^ADM- 

Equations (98) with (96) immediately give us the 0{t]) 
piece of Tadm ^ fixed angular velocity: = 

— {r]/2)dzi/dil, where we introduced the operator to 
denote a linear variation with respect to r] at fixed X. To 
obtain the shift at fixed energy, (5Tadm — 
write 




— TjP 

— P P^ADM 


dn 

'^^ADM y(n) mP 
jpP " P^ADM 
“P^ADM 




(99) 


where in the second line we used = 

— {d^l^ADM^ in the third line we 

applied dL^Y,M/'^^ADM ~ which is valid 

for any circular geodesic (omitting subdominant 
terms in 77 ). From Eqs. ( |9^ and ( |M| ) we find 
_ (^r]/2)[zi — Q{dzifdQ)], and substitut¬ 
ing this with the above result for arrive 

at the simple expression 


^^ADM- 


( 100 ) 


Note that in the analysis leading to Eq. (100) we have 


not assumed anything about the spin a of the central 
black hole, so the result should apply in general (and 
suggests an interesting new interpretation of Zi in terms 
of a shift in the critical value of the angular momentum). 
In the extremal case, D = l/2-|-0(e), so at leading order 
in e we obtain 


^^adm(^^) = -^^i(^^). (101) 
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where 


Zi{E) ■— limzi(f2(£;;e),e), 

e—)-0 


( 102 ) 


with the limit taken at fixed energy E. Here, for clar¬ 
ity, we have made explicit the functional dependence of 
zi and H on e, and have parametrized the circular orbits 
by their geodesic energy E, noting that the diffe rence be¬ 
tween E and subdominant in Eq. (1011. Indeed, 

in practice, Z\ (E) may be evaluated by considering a se¬ 
quence of circular geodesics of diminishing e (and a fixed 

E)- _ _ 

to ob¬ 


Equation (101) may now be used with Eq. 


tain the sought-for relation 


SL’^^^%E)=E^ + l-Zi{E), 


(103) 


which may then be used in place of (70) in both condi¬ 
tions (69) and (81). The relation (103) relieves us from 


the need to restrict attention to particles coming from 
infinity, which is why we have used in it the argument 
E in place of Eao. The energy may now t ake any v alue 
E > Eisco{= l/\/3), and the conditions (69) and (81) 


may be evaluated for all corresponding orbits. 

We may also use Eq. (101) directly in conjunction with 


Eq. (73), to write the censorship condition (in the generic 


case) in the remarkably simple form 


Z,{E)>-{E^-1). 


(104) 


Overspinning is averted if and only if this inequality holds 
for all E > Eisco- The evaluation of the condition (104) 


requires only redshift information on unstable circular 
orbits, evaluated at the extremal limit e —> 0 with E 
held fixed. 

It should be emphasized that the applicability of the 
theoretical framework underpinning Eq. (95) is yet to be 


rigorously established within our black-hole-perturbative 
context. There is, however, accumulating evidence to 
suggest it. The first-law framework has been tested ex¬ 
haustively against results from post-Newtonian theory 
in both the Schwarzschild [13] and Kerr m cases, and, 
more importantly for us, it has been shown to precisely 
reproduce certain rigorous GSF results in the strong-held 
regime, at least in the Schwarzschild case [43|. Further 
reassurance is provided, in the Kerr case, by the agree¬ 
ment between the hrst-law framework and the perturba¬ 
tive Hamiltonian on e I16| . Ideally, should be eval¬ 

uated via both Eq. (103) and the more rigorous GSF for¬ 


mula (70). Indeed, if nothing else, a demonstrated agree¬ 


ment between these two expressions would lend a strong 
support to the validity of both hrst-law and Hamiltonian 
frameworks in the strong-held regime. 


VI. ASSESSMENT OF NUMERICAL TASK 


A direct evaluation of the censorship condition (69) 


requires the function (5L“"®(Eoo), which involves the 


conservative piece of the GSF along critical equatorial 
geodesics on a Kerr background, in the extremal limit. 
To explore the case of hne-tuned orbits, via the condi¬ 
tion (81), would require, in addition, a calculation of the 
radiative huxes (either at inhnity or down the event hori¬ 
zon) for unstable circular equatorial geodesics, again in 
the extremal limit. Finally, to evaluate the overspinning 
condition in its alternative form (104) demands redshift 
data for these same unstable circular geodesics. To the 
best of our knowledge, these numerical tasks are beyond 
the capability of existing GSF codes—though, perhaps, 
not by much. We think that custom-built codes to pro¬ 
duce the necessary data could be developed through rel¬ 
atively mild adaptations of existing codes. In this section 
we review relevant numerical methods that are currently 
available, and discuss how they should be customized. 


There now exist two computational frameworks for 
strong-held GSF calculations in Kerr geometry. One, 
by Dolan and collaborators [IS], tackles the metric per¬ 
turbation equations in the time-domain (TD), via a di¬ 
rect numerical evolution of the hyperbolic set of lin¬ 
earized Einstein’s equations in the Lorenz gauge. A ju¬ 
diciously designed “puncture” scheme is applied to ex¬ 
tract the correct regular piece of the metric perturba¬ 
tion, from which the local GSF is calculated along the 
orbit. The second framework, due to Shah and collabo¬ 
rators gSlIlT], is based on Teukolsky’s perturbation for¬ 
malism: The relevant perturbation equations are decom¬ 
posed into Fourier-harmonic modes and tackled numeri¬ 
cally mode-by-mode, in the frequency domain (FD). The 
GSF is then reconstructed using a recently formulated 
mode-sum regularization procedure jUj. Both methods 
are underpinned by (the same) rigorous theory, and have 
been tested against each other and against results from 
post-Newtonian theory. All applications so far consid¬ 
ered the GSF on fixed geodesic orbits, neglecting back- 
reaction on the orbit. 


The above two methods are somewhat complementary 
with respect to the range of problems they can tackle. 
The FD method does best with bound-orbit configu¬ 
rations, where the perturbation field admits a discrete 
spectrum. It is not immediately clear how to apply 
the method to orbits that come from infinity. The TD 
method, on the other hand, can handle all types of orbits 
equally well. The special case of circular orbits can be 
tackled by both methods, but much more efficiently in 
the FD, thanks to the simple spectrum of radiation from 
such orbits. 


Let us consider first the more straightforward of the 
aforementioned numerical tasks: calculation of energy 
and angular-momentum fluxes from circular orbits (for 
the fine-tuning case). This is a standard calculation, 
most efficiently performed by solving Teukolsky’s equa¬ 
tion in the FD, either numerically [IS] or via the semi- 
analytical method of Mano, Suzuki and Takasugi [50] . 
Usually calculations focus on stable circular orbits as 
the source of radiation, but exactly the same techniques 
should be applicable without change to unstable orbits 
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below the ISCO; the stability properties of the orbit are 
immaterial here. The only potential complication comes 
from the need to evaluate the fluxes in the extremal limit. 
Care would need to be taken in correctly identifying an 
inner “wave zone” in which to impose the boundary con¬ 
ditions, for each finite value e <?; 1. It may prove con¬ 
venient to work with a rescaled radial coordinate [such 
as f := (r — i?oh)/e] in order to better resolve the near¬ 
horizon wave dynamics. However, we do not see any 
issues of principle to hinder such a calculation, and it 
could be based on any of the existing platforms, such as 
the one by Shah et al. m- 


We next turn to the calculation of 5L^^^(E) using the 
form (103). This requires an eva luation of the redshift 
function zi(H;e) [recall Eq. (102)] on unstable circular 
orbits, at the limit e —?> 0. Calculations of zi for sta¬ 
ble (or otherwise near-ISCO) geodesic orbits have been 
performed using both FD [III [47] and TD [16] methods, 
in a non-extremal Kerr geometry. In the Schwarzschild 
case, such a calculation was performed (in the FD) even 
for unstable orbits, reaching very near the light ring m- 
The challenge is to extend these calculations to the near¬ 
extremal Kerr case. The issues here are similar to the 
ones affecting fluxes. In the FD, one would need to care¬ 
fully set inner boundary conditions, and also carefully 
monitor the convergence of the multipole mode-sum, par¬ 
ticularly at large energy (lessons can be learned from the 
Schwarzschild analysis of ED). Some development, tests, 
and a careful error analysis would be required, but the 
problem seems to us perfectly tractable. 


If one is satisfied with the level of rigour provided by 
the first-law and Hamiltonian formulations, then no fur¬ 
ther calculations would be needed: The question of over¬ 
spinning, for both generic and fine-tuned orbits, can be 
answered based on numerical data pertaining to circular 
orbits only. However, to establish a full confidence in 
the results, a direct evaluation of via Eq. (70) 

would be required. Since this involves GSF data on un¬ 
bound orbits, a TD treatment is preferable. So far, TD 
calculations of the GSF in Kerr have been limited to cir¬ 
cular (and equatorial) orbits in a non-extremal geometry 
[52j . However, it should not be hard to implement orbits 
that arrive from infinity, including critical orbits in the 
extremal case, at least in principle. To achieve this, cer¬ 
tain technical details would need to be addressed. In the 
rest of this section we give an assessment of the challenges 
involved. 


First, to probe the extremal limit, we would need GSF 
data for a sequence of spacetimes approaching extremal- 
ity. In the current code [S2| , accuracy has been observed 
to degrade with increasing a. In the circular-orbit case 
the code can handle spins up to a = 0.7-0.8M, beyond 
which the loss of accuracy is rapid. The cause of the accu¬ 
racy degradation is not quite clear yet, but work towards 
a full diagnosis of the problem, and towards its resolu¬ 
tion, is now in progress. It appears likely that the code 
could be pushed to very high spins with only moderate 
effort. 


Second, one would need to implement plunging 
(geodesic, critical) orbits. In a TD framework this would 
entail only a minor coding effort; the basic code architec¬ 
ture should remain intact. Two of the technical details 
that would need addressing are (I) the handling of the 
auxiliary worldtube (see [45] for details), which would 
now be moving to track the radial motion of the parti¬ 
cle; and (2) the treatment of “junk radiation” (again, see 
[IS] for details), whose problematic effect is expected to 
be more pronounced than for circular orbits. A method 
for tackling the latter problem has been developed and 
successfully implemented in the electromagnetic case, by 
Zimmerman et al. [5]. It remains to test its performance 
in the gravitational case. 

Perhaps the most significant remaining problem is that 
of the m = I mode instability—an issue identified and 
thoroughly analyzed in [45] . In numerical experiments, 
this particular azimuthal mode of the Lorenz-gauge met¬ 
ric perturbation appeared to develop a linear instability 
at late time, which so far could not be cured. The seed 
of the instability appears to be a certain non-radiative, 
Lorenz-gauge mode, which is perfectly regular on each 
time slice. Various methods have been tried in attempt 
to filter that mode out of the numerical solutions, so far 
without much success. A simple filter can be applied in 
the circular-orbit case [IS], giving a satisfactory ad-hoc 
solution, but the method would not work for non-circular 
orbits. We are aware of at least two parallel efforts, by 
two groups, to obtain a more fundamental solution to the 
problem, and we remain optimistic that the issue will be 
resolved soon. 

The problem of m = 1 mode instability is entirely 
avoided within a third computational framework, now 
being developed by a collaboration involving one of the 
authors (LB). The new method combines the simplicity 
of the Teukolsky formulation with the utility and flexi¬ 
bility of the TD approach. It essentially involves a TD 
implementation of the Teukolsky equation, together with 
a scheme for constructing the GSF in a certain (non- 
Lorenz) gauge. When completed, the code will offer a 
most natural tool for performing the calculation required 
here. 


VII. SUMMARY AND CONCLUDING 
REMARKS 

We studied the scenario in which a massive (spinless) 
particle is dropped into a nearly-extremal Kerr black hole 
on an equatorial-plane trajectory. For this scenario, we 
presented a systematic analysis of the censorship condi¬ 
tion at first order in the self-force approximation. One of 
our main goals was to determine what GSF information, 
precisely, would be needed in order to provide a definitive 
answer to the question of whether an over-extremal black 
hole was a possible outcome (within the classical theory). 
We achieved this by formulating concrete, necessary and 
sufficient conditions for overspinning to be averted; these 





20 


are given explicitly in terms of GSF quantities, ready for 
numerical evaluation. 

Along the way, we have established several interesting 
results: 


• When the GSF is ignored, overspinning can be 
achieved in a certain open domain of the parameter 
space, mapped here precisely for the first time. 


Overspinning is possible (when the GSF is ignored) 
only with particles thrown in from infinity. For any 
value of the initial energy-at-infinity, overspinning 
can be achieved by choosing the particle’s rest mass 
and angular momentum from within certain open 
intervals, as prescribed in the last sentence of Sec. 
|n) below Eq. 


• In the full-GSF case, a sufficient and necessary 
censorship condition for “generic” orbits may be 
formulated on the one-parameter family of critical 
geodesics. That condition is sensitive only to the 
conservative piece of the GSF. 


• A more general condition may be formulated to en¬ 
compass also fine-tuned orbits (ones whose parame¬ 
ters are exponentially fine-tuned to produce an adi¬ 
abatic evolution along the global attractor). That 
condition involves also the radiative fluxes from un¬ 
stable circular geodesics. 


• The conservative GSF effect may be reformulated 
in terms of the “redshift” variable. This results in 
an alternative form of the censorship conditions, 
which involves only perturbative quantities (red- 
shift and radiative fluxes) calculated on unstable 
circular geodesics. 


Our main results are expressed in Eqs. (69), (81) and 


(104). Equation ( |69| [with (70)] is the censorship condi¬ 
tion for generic orbits. It required as input 5L“"®(Eoo), 
the conservative GSF correction to the critical value of 
the angular momentum at a fixed E^o ■ The condition is a 
sujficient one in the sense that its validity for all E^o > I 
would imply that censorship is protected in our scenario, 
at least for sufficiently small values of the particle’s rest 
mass (this last caveat is to remind us that our analysis 
and conditions are formulated within the first-order self¬ 
force approximation). Equation ( |69[ ) is also a necessary 
condition, in the sense that its violation for any value 
Eoo > 1 would mean a direct infringement of cosmic cen¬ 
sorship. 


Equation (81) [with (74) or (79)] is the more general 
censorship condition that covers also Hne-tuned orbits. 
Its evaluation requires, in addition to also the 

radiative quantities £+ and ^(qc)) 

sociated with the adiabatic evolution along the global 
attractor. Without fine-tuning, and vanish at 

relevant order, and (81) reduces to the generic condi¬ 
tion ( [^ . With fine-tuning, the condition (811 should be 
evaluated on the two-parameter space of initial and final 


energies, {Eoo,Ef}. It is both sufficient and necessary, 
in the same sense as (69). 

Finally, Eq. (104) is a reformulation of (69) in terms of 
the redshift variable Zi{E), calculated on unstable circu¬ 
lar geodesics in the extremal limit (taken with fixed E). 
The more general condition (81) may also be formulated 
in terms of the redshift, by substituting for using 

Eq. (1031. This alternative form is more readily amenable 
to numerical evaluation, because it requires only circular- 
orbit information. An additional advantage is that ( |104 ) 
is applicable to any value E^o > Eisco of the initial en¬ 
ergy, without the restriction that particles have to be sent 
from infinity. However, the redshift formulation relies on 
some layers of non-rigorous theory, so a direct evaluation 
of the conservative GSF effect, via Eq. (70), would be 
desirable as a check. 

We are unable to predict, in advance of a numerical 
calculation, whether our censorship conditions (691 or 
(1041 hold. We are not familiar with an argument to 
suggest even the signs of the terms and Zi ap¬ 

pearing in these conditions. However, if we assume that 
the generic-orbit conditions (69) [or (104)] are satisfied, 
then it is possible to conclude that censorship is pro¬ 
tected also for fine-tuned orbits [i.e., the inequality (81) 
holds], assuming only a certain plausible lower bound on 
the flux ratio TZ{E). Specifically, one can show, with the 
help of Eqs. (84) and (89), that (81) is satisfied for all 
Eaa > Ef > Eisco if > — 3 - This lower bound 

lies comfortably below the ISCO value TZ{Eisco) ~ —0.1 
estimated in Ref. [53) . and we expect unstable circular 
orbits to be less superradiant than the ISCO [i.e, have 
7?.(i?) > TZ{Eisco)] on account of their frequency be¬ 
ing larger than the ISCO frequency (recall also TZ > 0 
for E > ^). If a numerical calculation of Tl{E) con¬ 
firms this expectation, it would mean that fine-tuning 
disfavours overspinning. In that case, establishing the 
simple inequality (69) [or (104)] would suffice for ruling 
out the overspinning scenario. 
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Appendix A: ADM energy and angular momentum 


This appendix gives a detailed derivation of Eqs. (44) 


and ( |45[ ). We consider two relativistic point particles in 
flat space, representing the black hole-particle system at 
an infinite separation. We are given the two rest masses, 
M and /i, and the particle’s energy ^Eao(> /r) and angu¬ 
lar momentum in a reference frame attached to the 
mass M. The mass M has an intrinsic spin Ma perpen¬ 
dicular to the orbital plane, and the mass /r is spinless. 
Our goal is to obtain the system’s total energy and an¬ 
gular momentum in the center-of-mass (CoM) frame. 

First, we note that, in the limit of infinite separation, 
both black-hole frame and CoM frame are inertial, and 
they are related via a simple Lorentz boost. Let x°‘ = 
(t,x,y,z) be a Cartesian frame centered at M, so that 
the spin of M is aligned with the z direction, and the 
particle’s orbit lies in the x-y plane. Denote the four- 
momenta of fi and M in the black-hole frame by 

= 0,0,0). (Al) 

The magnitude of particle’s 3-momentum satisfies 


\pa^ = ipi)r+ipl)r=AEi-i). 


(A2) 


The CoM system x°‘ is related to x°‘ via a Lorentz boost 
A“^ in the x-y plane. Let and = 

^'^yP{M) denote the 4-momenta of p and M in a;“. The 
CoM condition, 


P(M) + P(/i) — 0, 


(A3) 


yields two nontrivial equations for the two boost pa¬ 
rameters (3^ = v^/c and = u^/c (where and 
yy are components of the boost velocity). One hnds 


/3* -l-pAoo) ^ for i = x,y. Thus, using (A2|, 


/? = [(/ 3-)2 + {pvfY/2 ^ 


1 + yEao 


(A4) 


where p = p/M. 

Now that we have at hand the boost A‘)g(/3), the rel¬ 
ativistic energies of p and M in the CoM system are 
obtained as pO^^ = and = A°pP^^M)- The 

sum of these two energies is the total, ADM energy of 
the system. A short calculation gives 

Aadm = P(^) +P(M) = M{1 + 2riEao + . (A5) 

This result is valid for any mass ratio p. For p <C 1, an 
expansion in p through 0(p) gives Eq. (44). 


To obtain the total angular momentum in the CoM 
system, we need first to relate the particle’s CoM posi¬ 
tion X(^) and black hole’s CoM position X(m) to their 
separation x in the black-hole frame. This is achieved by 
solving the CoM condition 


Hp)P\t,) + ^(m)pIm) = 


(A6) 


simultaneously with |x(,i)| + |x(m)| = |x|, where i* = 
A^l^x^ {i = X, y) is the separation in the CoM frame. The 
particle’s CoM orbital angular momentum is then given 
by L(fi) = - y(fi)Pf^y and similarly for the mass 

M. The total, ADM angular momentum of the system 
(with respect to the CoM) is the sum of L(^), L(m) and 
the spin angular momentum. A short calculation, using 
Eqs. (A2) and (A4) and the relation pLoo = xpj^^ — yp^^y 
gives 


Tadm = M {a + pLoo/Eadm) ■ 


(A7) 


This is valid for any p. Substituting for Eadm and ex¬ 
panding through 0{Tf) produces Eq. (45). 


Appendix B: Shift in the critical value of the 
angular momentum due to the full self-force 


In this appendix we derive Eq. (51), which describes 
the GSF-induced shift JLoo in the critical value of 
the angular-momentum-at-infinity for a given energy-at- 
infinity, at leading order in p and in e. In Section [ill D I| 
we obtained an expression for [Eq. ( |4^ ] under the 
(non-physical) simplifying assumption that the GSF has 
no dissipative piece. Here we restore dissipation and cal¬ 
culate the shift 5Lao caused by the full GSF. 

In our analysis we only require the leading order [0(p)] 
of 5Lao- As explained in the text, for that purpose it is 
sufficient to ignore fine-tuning and assume the changes 
in (specific) energy and angular momentum along the 
critical orbit are O(plnp) at most. 

Consider, then, a critical orbit parametrized by Eao{> 
1), and an arbitrary moment t = tq after the orbit had 
settled into quasi-circular motion, but before AE{tq) and 
AL(to) have accumulated 0(1) changes [specifically, as¬ 
sume AE{tq), AL{tq) = O(plnp) at most]. Evaluate 
the full-GSF radial equation of motion (35) at time tq. 


subject to the near-circularity condition dr/dr = 0{r]), 
substituting for E{to) and L{to) from Eq. (37), and re¬ 
placing Lao with Lo{Eao) + SLao{Eao)- At leading order 
in p and in e one obtains the solution 


SLaoiEao) = AlE(ro), 


(Bl) 


where AW is the GSF integral defined in Eq. (50). This 


result can only make sense if (i) the expression on the 
right-hand side is in fact independent of tq at leading or¬ 
der; and (ii) the quantity AIT(ro) = 2AE{to) — AL(to) 
remains of 0(p) even for tq large enough that the individ¬ 
ual terms AE{to) and AL(to) are already of O(plnp). 
We now argue that both conditions are satisfied. 

To make the argument, let us split AIT(ro) into an 
“approach” piece AIT(rc), and a “quasi-circular” piece 


AtT(rc,To) := AIT(ro) - AIT(tc) 


-p ^ f {2Et + E^) dr. 

J Tr. 


(B2) 
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Here the end-of-approach time Tc can be taken to be any 
moment Tc < tq after the orbit had settled into quasi¬ 
circular motion, in the sense that r(Tc) = 0{ri) at most. 
We will now show that 


AWItcTo) = 0{e)0{r)lnT]) 


(B3) 


(or smaller). 

Next, consider the contribution AW(^t)' 

fTa pE(to) 

= (e/v) / b{E)FtdT =-e / b{E)dE 
Jtq J E{Tc) 

E{to) 


at most, for any choice of tq and of Tc- Assuming 
e| Inryl ^ 1, this would mean that AW(to) is dominated 
by its approach piece AW(tc), so that (i) AW(ro) does 
not depend on tq at leading order, and (ii) AW(ro) ~ 
AW{tc) = 0{r]) as argued ab ove. 

To establish the scaling in (B3), use the orthogonality 


= - eV6A2 - 2 


E{r^) 


(B7) 


In the second equality we have used Ft = —/r dE/dr, and 
in the third we have substituted for b{E) and integrated 
explicitly. Since the energy difference E{tq) — E{tc) is at 
most of 0{r]\nr]), we conclude that 


relation u“T’q, = 0 to write Ft + flE^ = —(u^/u*)Fr., 
where B and the replacement u°‘ —>■ does 

not affect the expression at leading order in rj. We wish 
to evaluate this relation for Tc < t < tq, when the orbit 
is quasi-circular. At leading order in 77 , B e) is 

then equal to the angular velocity of an unstable circular 
geodesic of specific energy E = E{t). At fixed energy, the 
angular velocity admits the small-e expansion B ^ — 


AW(t) = 0{e)0{r]\nr]) 


(B 8 ) 


(or smaller). 


The combination of Eqs. (B 61 and (B 81 leads to the 


\b{E)e -H 0(e2) with b = 6E{6E^ - 2)-^/^ [see Eqs. (87) 


and ( 88 )]. Thus, omitting terms that are subdominant m 


scaling stated in Eq. (B3). The upshot is that the 
contribution to AIT(ro) from the quasi-circular part, 
AW(Tc,ro), is negligible compared to the contribution 
from the approach part, AW(tc) = 0 (e°) 0 (? 7 ) (assuming 
ejlnTyj 1) . In other words, the GSF integral AIT(to) 
in Eq. (iBlj) may be truncated at the end-of-approach 


the integrand in Eq. (^| is 2Ft + F^ — 2{u'^/u*)Fr+ time Tc, with the latter taken to be any instance after 

\beFtj,, or, equivalently. 


2Ft +F^ = -2{u^/u*)Fr - beFt. 


(B4) 


the orbit had settled into quasi-circular motion [but be¬ 
fore the specific energy has accumulated 0 ( 1 ) changes]. 
This establishes Eq. (51) in the main text. 


Let us denote the contributions to AIT(to,Tc) from the 
first and second terms on the right-hand side of (B4) 


by AW(r) and AW^t}y respectively. In what follows we 
consider each of the two contributions in turn. 

Start with AW(^r) ; given by 


fri-ro) 

AIT(r) = ( 2 / 77 ) / (Fr/u*) dr. 


Note r(ro) — r(rc) = 0(e), since both radii belong to un¬ 
stable circular geodesics. From Eq. (43) we recall that 
for |f| <?; 1 the radial component Fj. is dominated by 
its conservative piece, ~ For(r), which is approxi¬ 

mately constant within the integration domain and may 
therefore (as we only keep track of the leading term in 
e) be pulled out of the integral. A simple calculation 
gives l/u* oc e (at fixed A), and this factor can likewise 
be taken out of the integral. We thus obtain the scal¬ 
ing AWt^r) ~ where Fq^ is evaluated, e.g., at 

r = r{Tc), the end of approach. It remains to determine 
the e-scaling of Numerical evidence, to be presented 
elsewhere [36], suggests the scaling Fo^. ^ e“^. This is 
consistent with what one would obtain by assuming that 
the GSF components in a normalized coordinate basis 
are finite: For = Fof ( 3 ’’’')“^^^ ~ e“^ (or smaller), assum¬ 
ing the normalized component Fof is finite and noting 
^ e^. Assuming, therefore, that Fq^ does not diverge 
faster than ~ e“^, and recalling Fg^ oc rf as usual, we 
arrive at 


Appendix C: Radiation from transition to plunge 
and final plunge 

In this appendix we argue that the term in 

Eq. (64) is subdominant for 77 , e ^ 1 and may there- 


(B5) fore be dropped within our leading-order analysis. Recall 
W(t„d) = 2 f(t„d) -^^end)’ ^here and are the 

energy and angular momentum radiated out to infinity 
during the transition from adiabatic inspiral along the at¬ 
tractor to a final plunge into the black hole, and during 
the plunge itself. 

Gritical orbits may transit into plunge in one of two 
ways: If the orbit is perfectly fine-tuned, the transition 
will occur around the location of the ISGO, and will then 
be similar to the transition at the end of a physical adi¬ 
abatic inspiral (on a stable orbit) that has already been 
studied in detail [321 IMl IMl IM]- With any less than 
perfect fine-tuning, the particle will slide off the peak 
of the effective potential and into the black hole before 
the ISGO is reached (cf. Fig. |^. It is reasonable to ex¬ 
pect the former scenario (transition through the ISGO) 
to yield the maximal radiation output, because (i) or¬ 
bits linger much longer around the ISGO location, where 
the potential is very flat, than they do around the peak 
of the potential; and (ii) the remaining distance to the 
horizon is maximal when the transition is at the ISGO. 
Below we start by looking at the “worse case” scenario 
of transition through the ISGO, and argue that, even in 


(B 6 ) 


that case, is negligible in Eq. (64). We then also 


AIT(r) = 0 (e) 0 (? 7 ) 
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examine the more generic scenario of a transition from 
an unstable orbit, and, as expected, arrive at a similar 
conclusion. 


1. Plunge from the ISCO 

For a transition through the ISCO, we use results by 
Ori and Thorne [32j , who studied radiation from the tran¬ 
sition regime at the end of a quasi-circular inspiral. An 
inspection of their analysis reveals that the main results 
are insensitive to whether the particle arrives the transi¬ 
tion regime along stable or unstable orbits. In particular, 
their equation (3.26) for the deficits in (specific) energy 
and angular moment over the entire transition should 
hold in either case. We write it here in the form 

(AL)tra„s = -A(a)^7^/^ 

('^-^)trans ~ ^isco A(a)^ 7 "/^ (Cl) 

where A(a) is a certain (positive) function of the spin a 
only, given explicitly in [32) . These expressions hold for 
any a, at leading order in ^ 1. The factor 77 ^/^ arises, 
essentially, from the fact that the transiting orbit spends 
an amount of proper time oc 77 “^/® whirling around the 
ISCO location on a nearly circular orbit [cf. Eq. (3.20) of 

The function A{a) involves a certain a-dependent di¬ 
mensionless factor, denoted £ in [32] (not to be con¬ 
fused with our £^), which describes corrections to the 
leading-order quadrupole emission formula for circular 
orbits, and is to be determined numerically. Kesden [33] 
used numerical data by Hughes (from a code presented 
in m) to estimate, in the near-extremal case, £ oc 
a scaling previously suggested by Chrzanowski [56]. We 
have confirmed this scaling using much more accurate 
numerical results by M. van de Meent 133, to be pre¬ 
sented elsewhere [36]. We also note that the scaling 
£ oc follows simply from a regularity assumption, 
namely that dE/dr must remain bounded (and nonzero) 
even in the limit e —?► 0: Noting that the Ori-Thorne func¬ 
tion £ is defined with respect to coordinate time t (not 
proper time r), and that {dr /dt)isco oc we obtain 
£ oc (dE/dr)(dr/d<)isco oc Assuming this scaling, 

and expanding the remaining a-dependence of A{a) in e, 
we obtain, at leading order in e, 

(AL)trans = , 

(AE)t,ans = (C2) 

with some (positive) numerical coefficient oq whose ex¬ 
plicit value will not be needed here. 

Refs. [33] [35] discuss the reason for the non-physical 
divergence of (AE)trans and (AL)trans when e —>■ 0 is 
taken with a fixed 77 , but this will not concern us here. We 
are instead interested in the combination (AIT)trans := 
2(AE)trans-(AL)trans, which, noting Oisco = |-fO(e^/^) 
and keeping only the leading term, reads 

(AIT) 

trans — 6oe«/l^774/^ (C3) 


with some (positive) numerical coefficient bo- Assum¬ 
ing the usual balance between the local dissipative GSF 
and the flux of energy and angular momentum in grav¬ 
itational waves emitted during the transition, we have 

Wians = -?7(^W^) trans and thuS 

W+ans = -boe^/^S^'^- (C4) 


Now examine the magnitude of Wtrans compared to 
that of other terms in the censorship condition (64|. If 
e > 77 , we have which, for e ^ 1, is 

much smaller than the term in Eq. (64). If instead 
we have e < 77 , then which, for 77 <C 1, 

is much smaller than the 0{q^) terms in that equation. 
The conclusion is that the contribution to from 

the transition regime, W/Zans, always subdominant in 
Eq. (64) for 77 , e <C 1. 

It remains to assess the contribution to from the 

final plunge into the hole. In general (when the black hole 
is not near-extremal) one identifies a final stage, extend¬ 
ing smoothly from the transition regime, where radiation 
reaction is negligible and the orbit plunges into the black 
hole on a nearly geodesic trajectory [321|33]. The picture 
may change a little in the near-extremal case, because 
the radial velocity remains small, |f'| oc e <C 1, all the 
way to the horizon. This can mean that the conditions 
that define the transition regime never quite break down 
before the the horizon is reached. In other words, there 
is a possibility that the particle crosses the horizon while 
still in the transition regime]^ 

That possibility can be assessed using Eq. (3.20) of 
[33], according to which the radial extent (Ar)trans of 
the transition regime, in the near-extremal case, is oc 
AA^-rf-A with a coefficient of order unity. This should 
be compared with the radial distance from the ISCO to 
the horizon, Ar ~ (2e)^/^ [recall Eqs. (10) and (13)]. It 
follows that for e /77 smaller than a number of order unity, 
the transition regime extends to the hori zon. In such 
cases, the amplitude of Wtrans ^q. (C4) serves as an 
upper bound for the amplitude of , and it follows 

immediately that the entire term is negligible in 

Eq. ([el. 


Let us then consider the case where the transition ends 
before the horizon is reached, so that a plunge regime is 
identifiable. In the plunge regime, by definition, the mo¬ 
tion is very nearly geodesic, and the near-horizon analy¬ 
sis of Mino and Brink [51] should be applicable. Ref. [31] 
obtained (among other things) an analytic expression for 
the energy output from the final plunge, by analyzing so¬ 
lutions to Teukolsky’s equation in the near-horizon, low- 
frequency, quadrupole approximations. They find (in our 


This possibility was studied in some detail in Ref. m ; see, in 
particular, Fig. 5 therein, in which ‘5’ is equivalent to our e^. 
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notation' iTp 


^I^ung ^ e (nnit - reh)iE - ^hL)- 


(C5) 


at leading order in e, where ri„it is the radius at the 
start of the plunge, and Teh is the horizon’s radius (de¬ 
noted i?oh in the main text). For a plunge following 
a transition at the ISCO, rinit — reh = and 

E-QhL ^ 0{e), giving Thus 

£j)(ung i® strongly suppressed at small e. Since the mo¬ 
tion is nearly circular even during the plunge, we have 
^Jiung - and the radiated angu¬ 

lar momentum is similarly suppressed. The combination 
Wpiung := 2£f^ung - -Cpiung 1® evcn more strongly sup¬ 
pressed at small e, and clearly contributes negligibly to 

W(t„d)- 

In conclusion, we have found that, for a transition 
through the ISCO, W(tndl= CS W^ans 1® 

always negligible in Eq. (64) for 7 y,e ^ 1. That radiation 


from the transition to plunge should have a negligible 
effect on the conditions for overspinning was previously 
suggested by Kesden |33] and Harada and Kimura [55] . 


2. Plunge from an unstable circular orbit 


This scenario is rather different from—and much sim¬ 
pler than—a transition through the ISCO. As the orbit 
is perturbed away from unstable equilibrium, its subse¬ 
quent evolution is almost immediately controlled by the 
“geosedic” radial force (proportional to the derivative of 
the effective potential), and back-reaction corrections be¬ 
come negligible. Let us state this point more precisely. 
Suppose r = rend is the radius at which the particle leaves 
the attractor (for concreteness, this may be chosen as the 
radius of the last turning point along the attractor). Note 
that the radial acceleration due to the geodesic potential 
is oc (r — Tend)) with a coefficient of order unity [since 
the second derivative of the effective potential at rend is 


of 0(e°)]. Since rend — ?'eh = 0(e), the magnitude of the 
geodesic radial acceleration is of 0(e) throughout much 
of the plunge. This should be compared with the magni¬ 
tude of the radial self-acceleration, which is of 0 (? 7 e)j^ 
We may define the onset of plunge as the point where 
the geodesic acceleration takes over from the GSF in con¬ 
trolling the motion; this happens near a radius r = rping 
satisfying rping —rend = 0(r]e). For r < ?’ping, the motion 
is governed by the geodesic equation of motion ([^ to a 
good approximation. 

We wish to bound the magnitude of suffi- 

cien tly well to show that it contributes negligibly in Eq. 
(64). Let us first consider the contribution to 
from the pre-plunge orbital segment rping ^ r < rend- 
The proper-time interval along this segment is At ~ 
(rend - rpins)/r = 0{r]°e°) (at most), since f = 0{r]e) 
along the attractor, where the evolution is driven by ra¬ 
diation reaction. Hence, the change experienced by the 
specific energy and angular momentum along this seg¬ 
ment is AE ~ HAL = 0(r]). Recalling H = 1/2 -|- 0(e), 
this gives AW = 0(er]) and thus a contribution of 0(erf) 
to negligible compared to the 0(rf) terms in Eq. 

(@. 

Next consider the contribution to from the 

plunge segment reh < r < rping. Using the geodesic 
equation of motion (§, one finds that the proper-time 
interval along that segment is cx ln[(rend — ?'eh)/(?'end — 
J'ping)] In 77 . The corresponding change in specific en¬ 
ergy is AE = 0(77 In 77 ). Since the radial velocity remains 
small, r = 0(e), throughout the entire plunge, we have 
AL ~ H”^AL, giving AW = 0(eri\ogri) and a contribu¬ 
tion of 0(6772 log 77 ) to Once again, this is negli¬ 

gible compared to the 0(ri^) terms in Eq. (64), assuming 
only e <C I/| In 77 / 

We conclude that, whether the plunge from the attrac¬ 
tor occurs at the ISCO or earlier, the term in Eq. 

(64) is always sub-dominant and negligible for e, 77 <C 1. 
Within our approximation, the energy the particle car¬ 
ries with it as it crosses the horizon is the energy with 
which it has left the attractor. 
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